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ABSTRACT 


Nonlinear  forced  oscillations  of  slender  beams  are 
studied.  The  analysis  takes  into  account  both  the  nonlinear 
effects  arising  from  large  deflections  of  the  beam  and  those 
arising-  from  nonlinear  material  behavior.  A  hysteretic 
stress-strain  law  of  the  Davidenkov  type  is  used  in  the 
analysis.  Detailed  results  are  given  for  large  amplitude 
oscillations  of  beams  with  hinged  ends.  Theoretical  results 
for  a  simply- supported  beam  are  compared  with  experimental 
results . 
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SECTION  1 


INTRODUCTION 


The  study  of  nonlinear  vibration  of  elastic  beams  and  plates 
has  had  a  long  history.  The  most  widely  known  large  deflection 
theory  of  plates  was  established  by  von  Karman  in  1910  £lj  ,  An 
account  of  the  derivation  of  the  von  Karman's  equations  is  given 
by  Timoshenko  etal.,  .  Biot  |^3a,  3bJ  formulated  the  theory  of 
large  deflection  of  plates  directly  from  the  nonlinear  theory  of 
elasticity,  Novozhilov  used  a  similar  approach-  He  made  the 
distinction  between  the  case  of  strong  bending  and  the  case  of  inter¬ 
mediate  bending.  In  the  former  case,  no  restriction  is  placed  on 
the  magnitude  of  the  rotation  of  the  cross  section;  while  in  the 
latter  case,  which  is  the  same  as  von  Karman’s  theory,  the  angle 
of  rotation  is  restricted  to  be  small  as  compared  with  unity. 

Based  upon  von  Karman's  theory,  a  number  of  problems  of  rec- 


5  a 


and 


,  5bj  and  Way  |^6a|  , 

6bj  and  Stippes  and  Haas  rath 


tangular  plate  were  solved  by  Levy 
problems  of  circular  plates  by  Way 

,  The  dynamical  equations  corresponding  to  von  Karman's 
theory  were  studied  by  Herrmann  and  more  recently  by 
Tadjbakhsh  and  Saibel  .  Eringen  £l0a,  lObj  ,  on  the  other  hand, 
formulated  the  dynamical  equations  for  beams  and  later  also  for 
membranes  which  correspond  to  Novozhilov’s  strong  bending  case. 
The  nonlinear  free  vibrations  of  string  has  been  studied  by  Carrier 
£lla,  llbj  and  using  methods  similar  to  that  of  Carrier,  the  non¬ 
linear  free  vibration  ui  beams  having  immovable  hinged  ends  has 


^Numbers  in  brackets  refer  to  listed  references  begin¬ 
ning  on  page  87. 
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been  treated  by  Eringen  Oaj  *  Both  Carrier  and  Eringen  use 
a  perturbation  method.  The  same  perturbation  technique  was 
extended  in  the  study  of  free  vibration  of  rectangular  plates  with 
hinged  immovable  edges  by  Chu  and  Herrmann  ^izj  .  As  for  the 
forced  vibration  cases,  approximate  solutions  for  the  nonlinear 
response  of  rectangular  plates  were  obtained  by  Kirchman  and 
Greenspon  j^l  3j  with  the  use  of  the  static  load-deflection  relation¬ 
ship*  The  nonlinear  response  of  a  beam  with  immovable  hinged 
ends  due  to  an  excitation  sinusoidally  distributed  in  the  space 
variable  was  treated  by  Mettler  ^14aj  ,  By  assuming  the  resultant 
axial  force  to  be  independent  of  the  spatial  coordinate,  he  was  able 
to  obtain  a  Duffing  type  response  showing  a  hard  spring  effect.  A 
similar  treatment  was  also  given  by  Kauderer  £l&J.  The  bending 
equation  used  by  Mettler  and  Kauderer  were  the  same  one  as 
originated  by  Kirchhoff  j^lbj  which  are  essentially  the  same  as 
the  one  of  the  von  Karman’s  theory.  By  assuming  a  first  order 
approximation  solution,  Yamaki  also  obtained  a  hard  spring 

Duffing  type  response  for  rectangular  and  circular  plates.  The 
response  for  circular  plates  was  also  studied  by  Nowinski  j^isj  . 

In  considering  the  effect  on  bending  due  to  axial  force,  an  approxi¬ 
mate  solution  in  terms  of  elliptic  functions  was  obtained  by 
Woinowsky-Krieger  j^l9|  •  McDonald  ^Z()J  has  obtained  results  in 
terms  of  elliptic  functions  for  the  case  of  immovable  hinged  bars 
with  small  motion.  Experimentally,  the  nonlinear  hard  spring 
Duffing  type  response  were  shown  to  exist  by  Lee  ^Zlj  ,  Tobias  ^Zzj  , 


Lassiter  et  al 


•’  [23] 


and  Smith  et  al 


~  M  • 


2 


In  all  the  above  work  the  nonlinear  effects  arise  from 
large  deformations.  Another  source  of  nonlinearity  arises 
from  the  nonlinear  behavior  of  the  material  of  the  beam. 

Resonant  structural  vibrations  can  be  controlled  by  the 
effective  use  of  material  damping.  The  problems  of  practical 
interest  are  those  connected  with  the  vibration  of  metallic 
structures.  As  has  been  shown  by  Lazan  ^Z5j  ,  Lazan  and 
Goodman  |^26j  ,  and  also  others,  most  metallic  materials  under 
cyclic  loading  exhibit  a  relationship  between  stress  and  strain 
that  is  not  elastic  even  at  stresses  well  below  the  yield  point. 
The  stress -strain  relationship  during  loading  and  unloading  are 
different.  The  area  enclosed  by  the  loading  and  unloading 
branches  of  the  stress -strain  curve  {hysteresis  loop)  serves  as 
the  indication  of  amount  of  energy  dissipated.  This  area  of  the 
loop  is  independent  of  the  frequency  of  cyclic  strain.  Due  to  the 
above  reasons,  the  classical  viscoelastic  Voight  models  cannot 
adequately  describe  the  behavior  of  metallic  materials.  In  view 
of  this,  Mindlin  etal.,  ^2?j  proposed  a  semi-empirical  method. 

In  this  work  the  viscous  coefficients  are  taken  inversely  pro¬ 
portional  to  the  frequency  of  the  excitation,  so  as  to  compensate 
the  frequency  effect  when  the  time  derivative  of  strain  is  taken. 
Recently  Sethna  J^fjj  ,  in  dealing  with  transient  vibrations  of  a 
beam,  also  adopted  the  same  concept  but  to  a  nonlinear  visco¬ 
elastic  stress -strain  law.  A  completely  empirical  hysteretic 
stress -strain  law  has  been  proposed  by  Davidenkov  £29]  *  A  pair 
of  equations,  one  for  loading,  and  one  for  unloading  are  used. 


Pisarenko 


M- 


as  well  as  Panovko 


M- 


used  Davidenkov1  s 


stress -strain  law  to  study  the  small  vibrations  of  cantilever 
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beams  with  results  that  compare  very  well  with  experimental 
results. 

It  thus  appears  that  there  is  a  large  amount  of  literature 
concerned  with  nonlinear  effects  on  beam  vibrations  where  the 
nonlinear  effects  arise  from  large  deflections  of  the  beam.  In 
the  works  concerned  with  these  problems  the  material  is  always 
assumed  to  be  elastic.  Then  there  is  a  smaller  body  of  literature 
concerned  with  nonlinear  effects  that  arise  from  nonlinear 
behavior  of  the  material  of  the  beam  and  in  these  works  the 
deflections  of  the  beams  are  assumed  small. 

The  present  work  is  concerned  primarily  with  those 
problems  where  the  nonlinear  effects  of  large  deformations 
occur  along  with  the  nonlinear  effects  due  to  nonlinear  material 
behavior. 

In  surveying  the  literature  on  large  deflection  problems 
for  elastic  beams  it  appears  that  most  of  the  investigations  are 
concerned  with  slender  pinned  beams  with  immovable  ends  and 
a  hard  spring  Duffing  type  response  under  sinusoidal  excitation 
is  predicted.  Relatively  little  attention  is  given  to  problems  with 
other  end  conditions  and  in  particular  nothing  is  known  about  the 
large  oscillation  of  an  elastic  beam  that  is  truly  simply  supported, 
i.e.,  with  one  end  free  to  move  in  the  direction  of  the  axis  of  the 
undeflected  beam.  A  secondary  objective  of  this  work  is  to 
investigate  this  problem. 

With  these  two  objectives  in  mind  and  for  simplicity  of 
presentation  this  dissertation  is  divided  into  two  parts.  Part 
one  deals  with  the  elastic  problem  involving  large  deformations 
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and  part  two  deals  with  the  large  deformation  problem  with 
nonlinear  material  behavior.  In  this  latter  part  Davidenkov's 
law  is  used.  In  both  parts,  two  particular  cases  are  treated. 

One  is  the  immovable  hinged  ends  case,  and  the  other  is  the 
simply  supported  case*  In  both  cases,  steady  state  solutions 
are  sought.  The  motions  treated  belong  to  Novozhilov’s 
strong  bending  type,  and  the  equations  of  motion  are  similar 
to  those  of  Eringen’s.  Both  rotatory  inertia  and  shear  de¬ 
formation  effects  are  included  in  the  analysis.  The  excitations 
are  either  in  the  form  of  a  periodic  force  uniformly  distributed 
or  in  the  form  of  a  periodic  motion  exerted  at  the  supports  and 
the  nonlinear  response  of  the  beam  is  studied. 

The  mathematical  problems  are  nonlinear  boundary  value 
problems.  The  method  of  analysis  used  here  is  the  perturbation 
method. 

The  stability  of  the  nonlinear  steady  state  solutions  are 
also  investigated  in  an  appendiX^and  here  an  adaptation  of  the 
asymptotic  method  of  Krylov,  Bogoliubov,  and  Mitropolsky 


The  appendix  referred  to  here  is  the  one  in  a 
doctoral  thesis  by  the  first  author  at  the  University 
of  Minnesota,  with  the  same  title  as  that  of  this 
report . 
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SECTION  2 


STATEMENT  OF  THE  PROBLEM 

Consider  a  homogeneous  rod  of  uniform  cross  section, 
with  a  longitudinal  plane  of  symmetry.  Let  the  right-handed 
Cartesian  coordinate  system  Oxyz  be  chosen  as  shown  in 
figure  1  with  Oxz  as  the  plane  of  symmetry  and  Ox  passing 
along  the  median  line.  The  length  of  the  rod  is  f 


Fig.  1  -  Beam  in  the  undeformed  state 
2.1  DEFORMATION 

The  following  three  assumptions  will  be  made: 

{a)  The  rod  is  free  of  load  in  the  y-direction  and  we  will 
consider  all  quantities  independent  of  y,  and  that  there 
is  no  displacement  in  the  y-di  recti  on, 

(b)  Plane  cross  sections  remain  plane,  though  they  are 

not  necessarily  perpendicular  to  the  median  line 

* 

during  the  deformation, 

( c )  Relative  elongation  in  the  z-direction  is  neglected. 


'  Prescott  in  deriving  Timoshenko  beam  equation,  made  use  of 
the  same  assumption.  See  M. 


6 


Now  examine  a  segment  of  a  length  Ax  taken  from  the  beam 
as  shown  in  figure  2.  During  the  deformation,  if  the  cross  section 
which  was  originally  at  x  and  perpendicular  to  the  median  line 
makes  an  angle  0  with  Oz,  and  if  the  shear  angle  is  denoted  by 
yt  then  a,  the  angle  which  the  median  line  makes  with  x-axis,  is: 

a  =  0  +  y  (I) 


Fig.  2  -  Displacement  of  the  beam  element 

JLet  a  material  point  with  coordinates  (x,  z)  in  the  undeformed  state 
have  new  coordinates  (x1,  zr)  in  the  deformed  state.  Then  with 
u(x,  y)  and  w(x,  y)  as  displacements  of  the  median  line  in  the  x  and 
z  directions  respectively,  we  have 

x*  =  x  +  u  “  z  sin  0 

z1  =  w  +  z  cos  0  (2) 

The  length  of  a  line  element  (dx,  o)  in  the  deformed  state  is  given  by 
ds  =  (dxr^  +  dz'^)2 


+  (z  cos  0), 


l2+[ 


1  +  u,  -  (z  sin  0) 


jo1 


dx 


(3) 


which  is  obtained  by  using  (2)  and  by  taking  limit  as  Ax  approaches 
zero.  A  subscript  after  a  comma  will  mean  differentiation. 
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(4) 


The  relative  elongation  in  the  x-direction  is  defined  as 


Equations  (3)  and  (4)  yield 

(5) 

(6) 

where  £  is  the  relative  elongation  of  the  median  line. 


€  =  £  -  z  9, 


with 


e  =  f  w. 


X 

2  2P 

+  (1  +  u,  ) 

X  X 


-  1 


From  the  geometry  of  the  deflection  curve  as  shown  in  figure  3, 
we  have, 


tan  a 


w, 

x 


1  +  u, 

X 


sin  a 


w, 

X 


i  +  e 


cos  a 


1  +  u, 

x 


i  +  e 


(7) 


where  £  is  as  defined  in  (6) 
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2.2  STRESS -STRAIN  RELATIONS 


If  o"  is  the  normal  stress*  and  c r  is  the  shear  stress, 

XX  xz 

then  by  Hooke’s  law 


o-  =  E  f  =  E  £  -  Ez  0 r 

XX  x 

cr  =  G  y 

xz 


Now  let 


N 
M  = 
Q  ^ 


-/■ 

’I': 

-/ 


dA  =  ea£ 


xx 

zdA  =  -EI0, 

xx  x 


tr  dA  =  k'GAy 
xz 


(8) 


(9) 


In  (9),  A  is  the  cross-sectional  area,  I  is  the  area  moment  of 
inertia  of  cross  section,  k1  is  the  shear-deflection  coefficient 
which  is  a  modifying  factor  as  used  by  Timoshenko £3 sj  and 
others,  the  value  of  which  depends  on  the  shape  of  the  cross 
section.  For  instance,  for  a  rectangular  cross  section,  k’  =  0.833 


2.3  EQUATIONS  OF  MOTION 


A  free  body  diagram  of  an  element  of  the  beam  is  shown  in 
figure  4.  The  element  is  bounded  by  the  top  and  bottom  surfaces 
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of  the  beam  and  on  the  two  sides  is  bounded  by  cross  sections 

of  beam  that  are  plane  and  which  were  normal  to  the  center  line 

before  deformation.  The  shear  forces  (Q)  are  parallel  to  these 

plane  elements  but  the  normal  forces  (N)  are  along  the  center 

line  and  are  not  necessarily  normal  to  these  planes.  It  will  be 

assumed  that  the  Cartesian  axis  system,  Oxyz  {reference  frame) 

can  have  an  acceleration  - — —  in  the  z  direction.  Then  from 

dt2 

figure  4  the  equations  of  motion  are: 


(N  cos  a),  -  (Q  sin  9),  =  p  A  u,  (10) 

x  x  tt 

(N  sin  a),  +  (Q  cos  0),  +  K(x,  t)  =  pA  (w  +  f(t)),  (11) 

XX  tt 

M,  -  Q  (1  +£)  cos  (a  -  0)  =  -J0,  02) 

x  tt 

where  K  (x,  t)  is  the  vertical  load  during  vibration,  p  is  the  mass 
density,  J  is  the  mass  moment  of  inertia  per  unit  length. 

There  are  eight  unknowns  occurring  in  equations  (10),  (11), 
and  (12),  namely  N,  Q,  M,  a,  0,  u,  w,  and  £  .  But  in  addition  to 
(10),  (11),  and  (12),  we  have  one  equation,  (6),  another  one  from 
(7),  three  more  from  (9).  Thus,  there  are  enough  equations  to 
take  care  of  the  unknowns. 

We  will  proceed  to  reduce  the  above  eight  equations  to  three 
equations  and  express  them  in  terms  of  u,  w,  and  0.  Using 
equations  (7),  (6),  and  (9),  equation  (12)  gives 


Q 


J  0 


’tt 


El  0, 

xx 


w,  sin  0  +  (1  +  u,  )  cos  0 

X  X 

. 


From  equations  (1)  and  (9),  we  have 


a 


0  +  y 


0 


Q 

k’  AG 


(13) 


(14) 
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From  equations  (7)  and  (9),  we  have 


N  =  EAC  =  EA 


1  +  u,. 


(15) 


cos  a 

Substituting  equations  (7),  (13),  (15),  into  (10)  and  (11),  we  obtain 

(  1  +  U’x 

u  +  eaJ  (1  +  uljc)  -  ,  2  ,211]  V 

^  w,  +  (1  +  u,  )  *  j 

[  \  x  J 


P  A 


J  6,  -  El  0, 

tt  XX 


w,  +  (1  +  u,  )  cot  0  (  ’  X 

X  X 


=  0 


(16) 


W, 


p  A  w,  +  EA  /  w, 
tt  \  x 


J  0.  -  El  0, 

tt  XX 

w,  tan  0  +  ( 1  +  u,  ) 

X  X 


2  ,,  .  »2 
w,  +  ( 1  +  u,  ) 

X  X 


+  K(x,  t)  -  pA  f(t),tt  =  0  (17) 


Substituting  equations  (14)  and  (13)  into  the  first  of  equation  (7), 
we  obtain 


tan  (  0  + 


J  0,  -  El  6, 

tt  xx 


k'  AG  Tw,  sin  0  +  (1  +  u,  )  cos 


f  1+"’, 


(18) 


x  x 

Thus,  we  have  three  equations,  (16),  (17),  and  (18)  in  terms 
of  three  unknowns,  w,  u,  and  0,  Up  to  now,  we  have  not  put  any 
restriction  on  the  magnitude  of  displacement  (u,  w)  and  rotation  (0), 
As  long  as  the  strains  are  within  the  elastic  limit,  these  three 
equations  can  give  results  for  values  of  u,  w,  and  0  of  any  magnitude. 


2*4  COMPARISON  OF  THE  EQUATIONS  WITH 
OTHER  KNOWN  THEORIES 

Either  equations  (19),  (11),  and  (12)  or  equations  (16),  (17), 

[ 

and  (18)  can  be  reduced  to  equations  corresponding  to  more 
simplified  theories* 
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2*4-l  Eringen's  Equations  of  Motion 

By  neglecting  the  shear  deformation,  i.e.,  by  letting 
6  =  a,  and  by  dropping  the  external  loads,  equations  {10),  (11),  and 
(I2)  can  be  reduced  to  Eringen's  equations  fl  0a|  . 


2>4'2  Equations  According  to  von  Kar man's  Finite 
Deformation  Theory 

If  we  drop  shear  deformation,  i.e.,  put  0  =  a,  omit 
rotary  inertia  term,  i.e.,  put  J  =  0,  omit  the  term  f(t),  and  consider 
€  as  very  small  compared  with  unity,  i.e.,  put  1  +  £  <=  lt  from 
equations  (10),  (11),  (12),  and  (15),  we  obtain 

-  pA  u,^  +  EA  [(1  +  u,^)  -  cos  0  j  -  (Q  sin  0),  =0  (19) 

-  pA  w,tt  +  EA  [{1  +  u,x)  tan  0  -  sin  0  ]  +  (Q  cos  0),  +  K(x,  t)  =  0 

-  EI  e,xx  -  °  =  0  (21) 

The  second  equation  of  (7)  becomes 


sin  0  =  w, 


(22) 


We  further  assume  that  0  is  sufficiently  small  so  that  we 
can  neglect  all  quantities  which  have  the  same  order  of  magnitude 
as  03 . 


Accordingly,  (22)  becomes 

0  =  w, 

x 

Substituting  (23)  into  (21) 


(23) 


U  =  “  ”•***  (2D 

Substituting  (24),  (23)  into  (19)  and  (20),  and  retaining  only  first 
and  second  order  terms*,  we  obtain 


- - - - — - - - - — - 

In  deriving  (25),  the  term  (Q  sin  0),  is  dropped  and  in  derivating 
(26),  tan  0  is  approximated  by  (0  +  03/3|j  ,sin  0  by  (0-03/3j)  and 
(Q  cos  0).  -  ~ 


by  Q,, 


(2  0) 
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(25) 


pA  u,  -  EA  (u,  +  \  w,  ^),  -  0 

r  tt  x  xx 

3 

pA  w,  -  EA  (u,  w,  +  i  w,  ),  +EIw,  =  K  (x,  t)  {26} 

K  tt  x  x  x  x  xxxx 

Equations  (25)  and  (26)  are  equations  which  correspond  to 

von  Karman's  theory  of  plateJand  are  also  the  equations  used  by 

Kauderer  j^l 5 ] - 


2.4.3  Linear  Theory  and  the  Timoshenko's 
Beam  Equations 

First  linearizing  the  second  of  equation  (7),  we 


obtain 


a  =  w, 

x 


Using  equation  (14),  it  becomes 

Q 


e  + 


k'  AG 


=  w. 


or 


Q  =  k'  (w,  -  0)  AG 


Now  linearizing  equation  (15),. 

N  =  EA  u, 


(27) 


(28) 


x 


Finally  by  using  equations  (27)  and  (28),  and  by 
linearizing  equations  (10),  (11),  and  (12),  we  obtain 


-  pA  u,  +  EA  u,  —  0 
H  tt  xx 

-pA  w,  +  Q,  =0 
tt  x 

pi  9,  -  El  0,  -  Q  =  0 

r  tt  xx 

We  have  dropped  the  external  load  term  and  the  term  f(t)  in 
equation  (30),  and  we  have  let  J  =  pi  in  equation  (31).  Equation 
(29)  is  the  wave  equation  for  the  longitudinal  waves. 


(29) 

(30) 

(31) 
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Substituting  equation  (27)  in  equations  (31)  and  (30), 


we  obtain 


-pi  0.  +  El  0,  +  k'  (w,  -  0)  AG  =  0 

tt  XX  X 

pA  w,  -  k1  (w,  -  0),  AG  =  0 

tt  XX 


(32) 

(33) 


Now  eliminating  0  in  equation  (32)  by  using  equation  (33),  we  obtain 

(34) 


pA  w,  +  El  w, 

tt  xxxx 


(pl+^)w,xxttt^-w.tttt=  0 


Equation  (34),  or  equations  (32)  and  (33)  are  called  Timoshenko's 

jJ c 

beam  equations. 


See 


Reference 


p.  331*  equations  (f  )  (m)  or  equation  (129)- 
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SECTION  3 


SOLUTION  OF  THE  PROBLEM 


3.1  EQUATIONS  OF  MOTION  IN  DIMENSIONLESS  FORM 


Equations  (16),  (17),  and  (18)  will  be  put  into  dimensionless 
form.  First  the  following  notation  is  introduced 


n 


ix)  2 

n 


pAi4 

E  I 


,  n 


1. 2,  3,  . 


(35) 


where  wn  is  the  natural  frequency  of  flexural  vibration  of  the  nth 
normal  mode,  and  Pn  is  the  corresponding  fcharacteristic  number, 
both  for  the  linear  case*. 

We  shall  be  interested  in  the  steady  state  solutions  of  the 
problem  under  conditions  of  sinusoidal  excitation.  It  is  not  difficult 


to  show  that  the  interesting  results  occur  when  the  excitation 
frequency  is  in  the  neighborhood  of  one  of  the  linear  natural 
frequencies  of  flexural  vibrations  of  the  system.  It  will  there¬ 
fore  be  assumed  that  the  frequency  of  excitation  is  in  the 
neighborhood  of  the  n**1  linear  frequency  The  symbol  £2n 

will  be  used  to  make  the  time  dimensionless,  thus: 


r  =  Slnt  (36) 

Instead  of  excitation  by  an  external  load,  the  periodic 
motions  of  the  system  can  also  be  caused  by  a  sinusoidal 
motion  of  the  support.  Thus  we  let 

f(t)  =  w  cos  £2  t 

„  n 

or  2~ 

f(t)  =  £2^  w  cos  £2  t  (37) 


¥See  reference  [35|,  p.  325  (116). 
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Now  let  us  introduce  the  following  dimensionless  quantities 
u  w  x  w  . 

u  =  r.w  =  =  wo  =7  (38) 

where  f  is  the  length  of  the  beam. 

Substituting  (35),  (36),  (37),  and  (38)  into  (16)  and  (17),  and 
multiplying  both  equations  by  f /AE,  we  obtain  the  equations  in 
dimensionless  form: 


p  4  n  2 

— —  — — r—  u,  +  /(l  +  u,f)  - 
A (2  <o2  TT  '  5 


1  +  u. 


i 


n 


jw,£2  +  (1  +  u,^)2 


1  i 

2 


*e 


.  r^f  2  e,  -  e,„ 

_I _ 1  EI  n  tt  4£ 

Af2\  '"’Z  +  (1  +  "’£)cot0  f 


=  0 


(39) 


p  4  2 

Inn 

— 7  - o —  w»  +W,- 

A)2  u  2  TT  1  « 


w 


*1 


n 


w,£  +  (1  + 


W 


•e 


Jf2  2 

t  ,  TT  n  e*  "  e’tt 
I  y  (EI  n  tt  4bi 


P  4n  2 

Ki  I  Pn  n 


+  a*2  I  w»t  tan  9  +  (1  +  u  t)  [»t  AE  A»2 
Ai  |  |  ,4  5  Ai  u 


n 


^  w  cos  t  =  0 
2  o 

(40) 


Similarly,  equation  (18)  becomes 


(1  +  u,^)  tan 


0  +• 


2  EI 

Jf2  0,  -  0’££ 

n  tt  |2  41 


'G  A^ 


w,  sin  0  +  (1  +  u,^  )  cos  0 


=  w,^  (41) 
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3.2  PERTURBATION  SERIES 


To  solve  equations  (39),  (40),  and  (41),  we  shall  develop  a 
perturbation  procedure.  First  we  assume  the  following  series: 

2  3 

W  =  Wj  q  +  W^  T|  +  T|  +  .... 

2  ,  3 

U  =  U  q  +  U  q  +  u  q  +  .... 

1  u  ^ 


2  2  2 

12  =  +  A,  q  +  A  q  +  .... 

n  n  In  2n 

2  3 

0  =  0  q  +  9  q  +  0  q  +  .... 

1  u 


(42) 


where  q  is  a  small  dimensionless  parameter,  0  <  q  «  1  depending 

on  the  slenderness  of  the  beam  as  discussed  below.  The  third  of 

equations  (42)  is  written  in  a  manner  such  that  SI  stays  in  the 

th 

neighborhood  of  the  n  linear  natural  frequency  of  flexural  vibration. 

2 

In  equations  (39)  and  (40),  there  appears  a  quantity  I/A f  ,  If  we 
2 

let  I  =  Ar  ,  where  r  is  the  radius  of  gyration  of  the  cross  section, 

2  2  2 

then  1/M  =  r  }i  ,  which  is  the  inverse  of  the  square  of  the 

slenderness  ratio.  In  the  following,  we  shall  consider  only  slender 
beams.  Thus  we  let 


I 


=  cq 


(43) 


where  c  is  a  dimensionless  constant  of  zero  order  in  q.  We  also  let 


E 

k'  G 


eq 


(44) 


To  express  0  in  terms  of  w  and  u,  equation  (41)  can  be  used. 

By  using  series  expansions  and  by  using  equations  (42),  (43),  and  (44), 
both  sides  of  equation  (41)  are  expressed  in  power  series  in  q,  and 
then  by  equating  coefficients  of  like  powers  of  q,  we  obtain 
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ei  *  wi.5 


0„  =  'v  -  u,  w  .  +  ew,  A>> 

2  2,4  1,4  1.1  1.444 


0,  =  w,  e  “  u,  t(w->  t  -  ui  t  w«  t  +  ewi 

3  3,.|  1,|  2,|  1,4  1,4  1,665 


w 


LI 


w  u 

1,4  2,4 


4  J  2 

*  ec  p  w  —  +  ew  tcc-e(u  >w  )  +e  w  (45) 

n  i,4tt  pi  2,444  i.l  1,1  .44  1.44444 

From  equation  (6)  and  first  of  equation  (9): 


lr=e  *  [»./  +  a  +  s.?)2 


1 

2  -  1 


(46) 


Using  approximations,  we  have 


N  r-  /  ,  21  2 

■sr=e  •  v*K{*‘  -i.t  h 


21  3 


+b.4  +  wuw2.r*  ui.i  *i.e  I" + 


From  equation  (13), 


(47) 


0., 


4  2 

T  T  P  ° 

I  J  n  n 


Q 

EA 


2  44  2  pi 

Af  \  Ai  K 


0, 


TTl 


U> 


n 


{1  +  l  ros  0  +  w, .  sin 

5  ^ 


Using  series  expansions  and  equations  (42),  (43),  and  (45),  we  have 


ST'”  (-cwi.e«,+”  (2cui.«wi.e6e‘cw2.«4  +  cui.eee'"i.e 


J  2  „  4 

cew  +— ~  c  0  w  .  )  +  - 

1.44444  pi  Kn  i«4  t  t 


(48) 
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Finally,  using  the  series  expansions  of  sine,  cosine  and  (1  +  x)  2 
and  substituting  equations  (42),  (43),  (44),  and  (4  5)  into  equations 
(39)  and  (40),  we  obtain 
4 


P 


■CTl 


n 


to 


2  2 
~  (co  +  A  1J+  ...)  (u  q  +  U  r]  + 
2  n  In  1 ,  tt  2,tt 


n 


[  2  2  2  2 

(l+  ui,£  71  +  uz  £n  +  +  ui  c  n  +  u->  **1  +  i  w,  ’i 


i.e1  2,e 

2 


i.e 


u,  e%2  +  -)J.e  +  =n‘»,.ee{’'Il£i"  +  ->-s  ‘  0 


(49) 


n  ,  2  .  .  2  - 

"  CT1  (<*>  +  a,„t1+  <w,  +  w,  ^  +  ...  '  w  cos  t) 

£  n  In  I  ,tt  2,tt  o 

co 

n 

f,  ±  2  ,  W  ,  2,  22  22 

+  (Wj  ^  T1+  w9  CT!  +  ...)  (U,  £  T1+  U,  £T)  +IW,  t  TJ  -  U,  ,  T, 


2.e 


1.6  1  Z,£ 


...)  ,  +  CT|  (CTl2  P  4  w  ~  -  w 

J  6  [  n  1.5TT  pi  1, 


1.6 
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w. 


+  ui.iei  wi.e’1  +  "i,«  wi.«e,r  -  "'i.eeee*,r  +-,( 


1/.  2  j  2  2  2  2 

+  ...)  (1  -  u  n  -  u  q  -|w  ti  +u  r\  +  ...) 
1  >5  6>b  * >S  * ib 


Ki 

,6  +  EA  ' 


Equating  coefficients  of  like  powers  ofq,  we  get 
1  n* 

”  '  ui.M-° 

2  4 

T]  !  6  w  +  w  =0 

n  l.rr  14666 


(51) 

(52) 


:f By  using  the  boundary  conditions,  we  will  show  later  that  U|  :  0, 
Consequently,  for  brevity,  we  will  drop  all  Uj  terms  henceforth. 


(50) 
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2 

ti  : 


0 


(53) 


(u  t  +  i  w  ), 

2,1  1,1  | 


3  4  p  4 

*  :  Pn  W2,tt  +  W2,£See  +  .  2  W1„tt"7(W1,6  U2,£ 


1 


w 


n 


4  J  *V 

*i.t  ’•« +  ewu««« ' c  pn  wi,eiTT7T *  -ET 


+  W  (3  COS  T 
o  n 


(54) 


’■*'  - c  C  u2.tt  +  (u3,s  twuw2.e,'i  +  c(wuwi.«e),r 0  (55) 


p 

4  4  n 

■n  :  p  w,  +  w  +  - T"  (A  w->  +  w,  } 

1  Kn  3,tt  2  In  2,tt  2n  1  ,rr 


n 

3  2 


—  (u  w_  +  —  w,  w  +  w  u _  £ ), £  +  (ew  ttt 
c  2  j  £,  Zf  ^  2  1  ^  2,£  lf^  3  £  2,^^^ 


+  e  Wl,£4444  *  3  W*»4  '  W1,C  U2,4  "  PI  2  eC  Pn  Wl,gTT,,§6| 

4 

f  4  J  C  Pn  Aln  _J_ 

|C  W2,£tt  pi  2  W1,£tt  pi 


0> 


n 


+  w.  ...  (u,  *  +  2  wit^] 


1,444  v  2,e 


e  w 


In  4 

— —  6  w  cos  T 
2  rn  o 


(56) 


n 


*  where  K  =  K  r^,  and  w  =  w  are  substituted 


o  o 
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,KA 

The  external  load  term  I — —land  the  term  due  to  the  sinusoidal 

4  '  El  / 

motion  of  the  supports  (w^P^  cost)  are  made  to  appear  in  equation 
(54).  In  this  way,  we  are  starting  the  perturbation  with  the  solution 
of  equation  (52)  as  the  generating  solution.  Thus,  the  generating 
solution  is  a  solution  of  the  free  vibration  problem. 

The  above  equations  are  applicable  for  all  types  of  con¬ 
ventional  boundary  conditions.  Note  that  since  the  axial 
displacement  u  is  also  considered  in  the  present  theory,  a 
distinction  must  be  made  between  the  conventional  simply  supported 
case  and  the  both  ends  hinged  case.  In  the  former  case,  one  end 
of  the  beam  is  allowed  to  move  axially;  while  in  the  latter 
case,  both  ends  are  immovable.  These  two  cases  will  be  worked 
out  in  detail.  The  technique  used  for  these  two  cases  can  also 
be  applied  to  other  cases,  although  the  analysis  for  other  cases 

is  much  more  complicated. 

3.3  BOTH  ENDS  HINGED  CASE 
The  boundary  conditions  are 

At  £  =  0;  1,  w  -  0,  u  =  0  and  M  =  0  (57) 

Since  M  =  -E  I  0  we  have  at  £  =  0;  1,  6  =  0 

,x  »A 

If  the  series  expansions  for  w,  u  and  0  of  equation  (42)  are  used, 
then  conditions  (57)  will  be  satisfied  if 

at  |  =  0,  1  w.  =  0 

u  =0  i  =  1,  2,  3  (58) 


*For  instance,  fixed  ends,  free  ends,  etc. 
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From  {51}  and  { 52 ) 


U\,ii  =  ° 


(59) 


4 

P  w  +  w  =  0 

n  1,tt  l, Hit 


(60) 


Equation  {59),  together  with  the  boundary  conditions  that  at  £  =  0,  1, 
uj  =  gives  the  following  solution 

uj  =  0  (61) 

Solving  equation  (60)  under  the  boundary  conditions  (58),  we  have 


n 


=  n  7T 


CD 


wi  = 


z 

m  =  1 


,  m 

a  sin  mirt  cos  -  t 

m  Z 

n 


(62) 


where  n  is  any  positive  integer. 
Since 


61  Wl.£ 


GO 


e' =  E 


.  m 

a  m  it  cos  rmrE  cos  - r 

m  2 

n 


(63) 


m  =  1 

Using  the  solution  (62)  for  wj,  and  taking  into  consideration  of  the 
boundary  conditions  that  at  £  =  o,l ,  u2  =  o,  equation  (53)  yields  the 
following  solution^ 

oo 

2 


1  2  -  I  V 

U2,£  2  Wl,£  "  4  ^ 


m  =  1 


2  2  2  2  m 

a  m  tt  ■  cos  - t  (64) 

m  2 

n 
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Equations  (61)  and  (62)  give  the  solution  in  the  first  approxi¬ 
mation.  This  solution  is  the  same  as  that  for  the  free  vibration 
of  the  linearized  problem.  In  order  to  proceed  to  the  next 
approximation,  it  is  advisable  to  take 

w  -  a  sin  nir£  cost,  n  =  1,  2,  3,  ...  (65) 

1  n 

r  =  Q  t 
n 


i.e.,  in  equation  (62)  we  let 

a  =  0,  m  i  n 
m 


a  =  a  ,  m  =  n 
m  n 


and  consequently  equation  (64)  becomes 


1  2  1  ,  .2  2 

+  —  w  =  —  (a  n  it)  cos  t 

Z  4  n 


(66) 


A  justification  for  taking  in  the  above  form  is  given  in 
Appendix  A.  We  turn  to  equation  (54) 


(n  tt)^  A 


(nir)  w  +  w  +  ' 

2,tt  2,£4£g 


to 


In  1  ,  13 

2-'V1,tt-c(U2.6W1.C+;W1,£  '■  | 


n 


.  t  K  f  .  . 

,  .4  J  o,44 

+  ew  -  c  ( n  ir )  w  .  —  =  -  +  w  n  tt  c  os  t  \  o  f } 

!,££££££  i,££tt  pi  ei 


Express  the  loading  as  a  Fourier's  sine  series 


3  no 

K  f  4 

K  =  — —  +  w  n  ir  cost  ”  f  \  B _ sin  mir4  jcosSl^t  (6B) 


El  o 


I 

m  =  1 


m 


n 


Refer  to  note  on  page  5. 
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Using  equations  (65),  (66)  and  (68),  equation  (67)  becomes 


GO 


^  »2.tt  + 


I 

m  -  1 


B  3in  rmrf;  cos  r 
m  ^ 


,  A  In  1  3  4  4, 

-  (nir)  a  - r  sin  nirg  cos  t  +  7—  a  n  w  sin  httE 

n  Z  1 6  c  n 

*1 

(3  cost  +  cos  3  t  )  -  (e  +  c  )  (n  tt)^  a  sin  mr£  cost  “  0  (69) 

pi  n 

In  the  right  side  of  (69),  the  terms  with  the  factor  sin  nrr| 

cost  are  the  secular  terms,  so  called  because  they  would  produce 

the  solution  for  Mw  11  of  the  form  tA  sinx  sin  n^,  which  then 

L  n 

would  destroy  the  periodicity  of  the  solution.  Just  as  in  the  theory 
of  oscillation  for  nonlinear  systems  ^3oJ,in  the  process  of 
eliminating  these  secular  terms,  we  will  be  able  to  find  the 
response  relation  of  the  system*  Following  Pisarenko  ^3oJ  we  let 
equation  (69)  be  multiplied  through  by  sin  mr£  cos  t,  and  then 
integrate  the  resulting  equation  throughout  the  whole  length  of  the 
beam  from  0  to  1 )  and  over  one  cycle  (r  from  0  to  Z tt ) . 

Note  that  by  using  integration  by  parts  and  by  applying  the 
boundary  condition  of  (58),  we  have  ' 

2tt  1  4  1 

//  (n  it)  w  +  w  \  sin  mr£  cos  t  d  £  d  t  =  0  (70) 

J  [  2,tt 


Taking  into  the  consideration  equation  (7  0),  we  have 
See  Reference  [30]  Chapter  IX,  Section  10,  p.  54. 
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2tt  1 


^  oo 


L,  TT  I  T - 1 

/  /  I 


A  A, 

_  ,  .  ,  .4  In 

B  sm  rmr§  cos  t  +  (mr)  a  - -  sin  nuf  cost 

m  n  2 


0  m  =  1 


co 


n 


I  3  4 

77  a  (n  tt)  sin  ntr£  (3  cos  t  +  cos  3  r) 
Ibc  n  * 


+  (e  +  c  —  )  (mr)^  a  sin  nir£  cos  r 
pi  n 


sin  nir £  cos  t  *  0  (71) 


Carrying  out  the  integration,  we  obtain 


4  A1 
B  +  (ntr)  a 
n 


_  7  TT  a  3  (*n)4  +  (e  +  c  ~){n  it)6  a  =  0 

n  c  lb  C  n  pi  n 

r 

n 


or 


co 


J_n 

2  16  c  n 


3  2  /  J  ,  .2 

a  -  (e  +  c  -7)  (nir) 

pi 


B 


n 


n 


/  \4 
(n  n)  a 


n 


(72) 


The  response  equation  can  thus  be  written  as 
2 


n  a 

n  .In 

2  “  +  2  T1"  2  ( pi  '  k’G 

cj  co  Ajt  1 

n  n 


1  I  J  ,  E  |  2 

+  7".  -  1  (n  it) 


3  At  ,  ,2 

TT-  (a  q) 

161  n 


B  q 

n 


2  q 


(n  it)4  qa  - 
n 


(73) 


As  might  be  expected,  the  linear  correction  terms  are  due  to  the 

J$£ 

rotatory  inertia  and  the  shear  deformation,  while  the  nonlinear 


*It  agrees  with  equation  (140)  of  reference  [3&j  ,  p.  335 
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3  AJ  2 

term  ■ —  (-r^a  )  is  due  to  the  axial  force,  the  magnitude 

161  ”  ^  1  Af2 

of  which  is  proportional  to  the  slenderness  constant  —  =  — —  = 

f2.  nc  i 

—  .  When  the  term  containing  the  external  force  (Bn)  is  omitted 

in  equation  (73),  the  resulting  response  equation  represents  the 

3  2 

stem  of  the  response  curve.  Since  the  term  ,  .  (a  t|)  is 

161  n 

proportional  to  the  square  of  the  amplitude  and  is  always 
positive,  the  stem  of  the  response  curve  bends  toward  the 


right  (figure  10).  This  is  so  called  hard  spring  effect.  In  fact, 

the  pattern  of  equation  (73)  is  similar  to  that  of  the  Buffing's 

* 

equation  of  nonlinear  oscillation  theory. 


Another  physical  phenomenon  associated  with  nonlinear 


spring  effect  is  the  jump  phenomenon.  Here  jump  refers  to  the 
phenomenon  of  a  sudden  and  abrupt  increase  or  decrease  of  the 
amplitude  when  the  excitation  frequency  changes  only  slightly. 

In  the  present  case,  it  happens  in  the  region  of  dfbc  of 
figure  12.  In  figure  12,  as  excitation  frequency  increases,  the 
amplitude  increases  from  point  a  up  to  point  b.  Further  increase 
of  the  frequency  will  cause  a  drop  in  the  amplitude  to  point  c.  On 
the  other  hand,  if  the  frequency  is  gradually  reduced  from  point  c 
toward  the  left,  the  amplitude  will  go  along  the  right  branch  of  the 
response  curve  up  to  point  f.  Further  decrease  of  the  frequency 
will  cause  a  sudden  jump  of  the  amplitude  from  point  f  to  point  d. 
The  portion  f  e  b  of  the  right  branch  of  the  response  curve  can  be 
shown  to  be  unstable.  (See  Appendix  B). 


* 

icicic 


Reference  p.  85  equation  (25)  and  p.  88  figure  23 

Reference  ^ 3 zj  p.  245  and  figure  80. 

Refer  to  note  on  page  5. 
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Let  us  examine  the  case  when  there  is  no  external  load  and 
the  beam  is  excited  by  the  sinusoidal  motion  of  the  supports  only,  Le., 


K  =  f f t ) ,  =  w  (nit)  cos  Q  t 

tt  o  n 


Expressing  w  in  terms  of  a  sine  series,  we  have 


4  3“ 

K  =  4n  tr  w  cos  Q  t 
o  n 


oo 

I 


sin  mrr£ 

m 


(74) 


m  =  1 ,  3 

Using  equation  (74),  the  response  equation  (73)  becomes 


n 


=  1 


U) 


A* 


n 


for  odd  n,  and 
2 


n 


=  1 


to 


n 


Ai 


:[  J  e  1 

[pi  +  k'G 

|  +  ■ 

_ 

_  ,  2-i 

„  2  4w  n 

~  (^a  )2-  - —  (75) 

2  n  nirna 
4h  nJ 


<->2(^-+ i^)+ Jz(ipn>Z  (76> 


for  even  n. 

The  difference  of  the  response  for  odd  n  or  even  n  in  the  linear 
case  will  be  discussed  in  Appendix  C.  To  solve  for  u^,  using 
equations  (65)  and  (66),  and  the  boundary  condition  that  at  £  =  0, 
u^  =  0,  thus  we  have 


1  2  2  2  2  .  2 
u  =  —  a  n  tr  cos  r(l  -  2  cos  nir^) 

Z.i  4  n 


*Use  is  made  of  the  Fourier  sine  series  formula,  valid  for 
f(£)  =_w  ,  -1  <  £  <  0,  and  f(£)  =  wQ,  0  <  £  <  1  (here  only  the  right 
half,  0  <  £  <  1  is  needed).  For  reference,  see  Churchill  "Fourier 
Series  and  Boundary  Value  Problems,"  McGraw-Hill,  1st  ed. 
p.  64  Prob.  3. 

Refer  to  note  on  page  5. 
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By  integration 


1  2 

u  =  — -  a  n n  sin  2nir|  (1  +  cos  2t) 

Z  16  n 


(77) 


In  view  of  equation  (72),  equation  (69)  becomes 


4  4 

nrrw  +  w 

2.TT  2,U£! 


OO 


L 


B  sin  mirti  cos  t 
m 


m=  1  ,m/n 


1  3  4  4, 

— —  a  n  tt  sin  ntrt  cos  3 r 
16c  n 


Solution  of  which  is 
on 


W2  = 


L 


m=  1  ,m^n 


B  sin  rmr£  cos  t  a 
m  n 

— — 4~1  +  TiI7  5m  n,r|  cos  3t 

(m  "  n  )  tt 


(78) 


Thus  u  and  w  can  be  written  as 


U  =  Uj  X]  +  U^  T) 


(79) 


-1  ,  2 

=  ~  (an)  n  tt  sin  Zn?Tt  (cos  ZQ  t  +  1 ) 
Ion'  n 


w  -  w  T[  +  w  Tj  =  (rja  )  sin  mrf:  cos  t 
o  1  Z  n  n 


oo 


z 


Ol  Bm)  sin  mir£  _  „  t 

- — - - - - -  cos  O  t 

4  4  4  n 

(m  -  n  )  tt 


m=  1  #m^n 
(qa  )3  A iZ 

n  ,  .  _  ^ 

+  - -  -  sin  nTTc  cos  3Q  t 

1281  *  n 


(80) 
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Once  ft  is  known,  the  amplitude  {r|a  )  can  be  determined  by 


n 


equation  (73),  and  the  solution  of  u  and  w  can  be  obtained  by 
equations  (79)  and  (80).  Equation  (79)  shows  that  period  of 
u  is  only  half  of  that  of  the  external  excitation  or  of  w. 

It  also  indicates  that  when  (1  +  cos  2ft  t)  =  0,  u  =  0,  and  thus 


n 

u  =  0,  when  w  =  0,  and  an  examination  of  equation  (79)  shows 

that  the  nodal  points  (u  =  0)  are  at  £  =  “■  which  is  the  middle 

12  ^ 

of  the  beam  and  at  £  =  —  Since  sin  Zniti  is  anti- 

2n  2n 

symmetrical  about  the  point  £  =  — ,  the  beam  thus  is  always 
stretching  out  or  compressing  in  a  symmetric  way  about  the 
middle  section. 


3.4  SIMPLY  SUPPORTED  CASE 

Assuming  that  the  left  end  is  hinged  and  the  right  end  is 
roller-supported,  the  boundary  conditions  are 


at 

£  = 

o, 

u  =  0, 

it 

o 

M  = 

0  (i.e.,  0, 

and  at 

£  = 

i, 

It 

O 

M  =  0 

CD 

=  0),  and 

£  F^  =  N  cos  a  -  Q  sin  0  =  0  (i.e.,  no  resultant  axial  force  at 
£  =  1 ).  If  series  expansions  are  used,  we  have 


at  |  =  0,  u.  =  0,  w  =  0,  0  =  0 

l  i  l,| 

i  =  1, 2,  3  ... 

and  at  4  =  w.  =  0,  0  .  =  0 

1  i.£ 

i  =  1, 2,  3  ... 


(81) 


(82) 
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Substituting  for  N,  Q  from  equations  (47)  and  (48),  we  have 


£  Fx  =  (N  cos  a  -  Q  sin  6)  =  EA  [qUj  ^  ^  (u^  +  j  w[ 


2) 


2.  4 


+”  (”3.e+wj,5w2.e+cwi.4'vi.ieru1.ewi.e  ,+“  - 

Thus  we  have,  at  £  =  1 , 


“1,5  ■  0 

“z,5  *  I  wi,52  =  0 


"3,e  +  wi,5w24  +  cwi.5wi,555  'ui.5wi,5  =°  <83) 

Equations  (81),  (82)^  and  (83)  are  the  boundary  conditions  for  the 
simple  supported  case. 

To  obtain  the  solution,  we  start  with  equations  (51),  (52),  and 

(53) 

Uue =  0 


P»  W1.TT  +  WlltK£  -  ° 

(“2,5+I  wu\*° 


Considering  also  the  boundary  conditions  (81),  (82),  and  (8  3),  we 
obtain 

u  —  0 

1 

w,  =*  a  sin  ntr£  cos  r 
In 

0.  =  a  n  tr  cos  nir£  cos  t 

In 

P^  =  nir,  t  =  a  t,  n  =  any  fixed  integer  (84) 
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Here  again  instead  of  a  series,  we  choose  to  be  a  one 

term  solution.  The  reason  of  this  choice  is  exactly  the  same  as 

that  in  the  case  of  immovable  hinged  ends* 

Solution  of  equation  (53)  together  with  the  boundary  condition 
1  2 

that  at  £  =  1,  u  +  —  w  t  =0  (equation  (83))  is 

2  j  £  2  A  t  s 


u  +  i  w  =  0  (for  all  r  and  4) 

2 ,  g  2  1 1  g 


(85) 


Thus 


and 


1  2 

uz,e  =  "  2  Wl,£ 

=  -  —  (nir  a  cos  nTr£  cos  r)2 
2  n 


12  1 

u  =  -  —  (a  nir)  {£  +  — -  sin  imr^)  (1  +  cos  2  t)  (86) 

2  8  n  2  it  n 


1  n 

To  obtain  — —  ,  we  again  use  equation  (54),  But  this  time 
co 

n 

we  have 

1  3 

(u  w  +  “  iv  )  =  0  (equation  (8  5)  ) 

Mb  ^  Mb 


Thus  instead  of  (73),  we  now  have  the  response  equation  up  to  the 
order  t/  as 


R 


2 


n 


2 


u 

n 


+ 


E 

k'G 


2  2 
n  tt 


q2B 


n 


4  4 

n  tt  rja 


(87) 


For  excitation  by  support  motion  only,  and  for  R^  =  uy 
(87)  becomes 


1 


(88) 


If  the  nonlinear  (a  ^)  term  is  dropped  and  for  n  =  1,  equation 
(75)  is  reduced  to  equation  (88). 


31 


As  is  remarked  in  Appendix  C,  this  formula  is  no  different 

from  the  response  relation  of  the  linear  case,  and  thus  we  conclude 

that  the  effect  of  large  amplitude  on  the  frequency  will  be  at  most 
2 

of  the  order  r|  . 

We  now  proceed  to  find  the  quantities  w_,  u_,  and  A  from 

2  3  2n 

equations  (54),  (55),  and  (56). 


For  the  sake  of  simplicity,  let  us  consider  the  case  when  n  =  1 
and  the  loading: 

oo 

—  2  2  V'  i 

K  =  it  w  cos  t  =  4ir  w  cos  t  ?  —  sinmnE  (89) 

°  °  _  i  i  m 

ITl  1  f  3  f 

Equation  (54)  becomes 


oo 


4  4.  ,  3 

ir  w  +  w  =  4tt  w  cos 

2 ,tt  2,£US  o  _  “  m 

m  1 1  i ,  ,< 


^  sin  rmr£ 


4  4 1  ,  .  6  . 

+  tt  — '-r  a  sin  it  %  cos  t  +  tt  (e  +  c)  a  sin  tt  £  cos  t 

Lt  L  1 


(90) 


But  equation  (87)  gives 


A  3  4  1 1  6  , 

4tt  w  +  tt  a  — —  +  tt  (e  +  c)  a  =  0 

o  1  Z  I 


U) 


1 


Thus  equation  (90)  becomes 


oo 


4  3 

TT  w  +  w  .  .  .  —  4tt  w  cos 

2,tt  2,££££  O  “  k  m 

m  j  ^  5  f  •  •  ■ 


y~  sin  m 


(91) 


Refer  to  note  on  page  5. 
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The  solution  of  (91)  is 


4  w 

< 

w  =  - 

2  TT 


00 


Esin  mirf 
- 2- 

4 

m  -  3 ,  5,...  m  (m  -  1 ) 


(92) 


To  find  we  use  equation  (55) 

u3.«  =  cPn4  U2.rr  •  (w2,e  Wl,|  +  CWl.e  wl.«t,’S 

Thus 


U3.S=CPn *  J  U2,TTd?-(w2.SWl.S  +  CWl,4Wl.eS?)+kl(T) 


(93) 


After  substitution,  it  becomes 


u 


3,£ 


6 

tt  2  /  £  cos  2tt£  . 

Fc  ai  ( t  •  |coa  2T  - a 


OQ 


.  I 


4  cos  m-rr^  cos  tt£ 


m  -  3,  5t... 


m4  -  1 


2  2  4  2  .  2 

COS  T  +  ca^  IT  COS  TT  £  COS  T  +  k,  (t) 


1 


By  further  integration,  we  have 


6  ,  i .  3 

u  =  i-  ca  2  i-.-SiHilS 

3  2  1  U  J  | 


OD 


cos  2r -  a 


i 


4  w  cos 
o 


1  —  4 

m-3,  5,**.  2  rr  { m  -  1) 


fsin  (m  -  1)  tt£  ,  sin(m+l)Tr£  .  2  4  2  /  £  sin2nfl 

“  ,  +  - T"i -  +  ca,  TT  COS  T  p-  +  - 2 

L  m  -  1  m+lj  1  |2  4  -TT  / 


+  (t)  +  k2(r> 
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Since  at  £  =  0,  =  0  and  at  £  =  1, 


it  yields 


and 


u3.i  +  w2.«  wi.t  +  cwuwi.e«'0. 


k2(T)  =  0 


VT)=  -  Vcii2  i  •  71  008  2 T 

1  4  ir  ' 


Thus; 


7T  2 

U3=TCal 


1  -  2  ir^  .  4^.  s  in  2  ir  4 

2  * +  6  3 


-  4  it 


8tt 


00 


cos  2  r  -  2  a 


.  I 

m  =  3  f  5 , , 


w  cos  r 
_ o _ 

m4  1 
m  -  1 


sin (m  -  1 )  -it  4  sin  (m  +  1 )  it 4 
m  -  1  m  +  1 


2  4  2  /  4  s  in  2  it  4 

+  C  a.  IT  COS  T  I  _  +  — - - *■ 


4  IT 


(94) 


A  w 

21  1 

To  obtain - r  ,  we  multiply  equation  (56)  by  —  d|d  r  and  integrate 

2  a  * 

ui  1 


as  before: 


r  i  p 


w . 


,2  "1,§ 

3,Tt  3,4eee  [  ueue  3  “2, 4-1.4’'' "2.444 


u  ,  w ,  +e  w 


~  j  4 

2  c  e  — 7  tt  w 


TT  TT 

4-  — -  A^,w,  +  —7-  A.,  w„ 


pi  1,£tt  /  2  Z1  1*tt  2  11  2)TT 


Cd 


1 


tl). 


1  1  3  2  , 

—  w  .  u^  +w  u  +  —  w  w 

C  I  J,4  3,4  2,4  2,4  2  1,4  2,41,4 


b> 


11  4 

—  IT  W  COS  T 

2  o 


1 


Jc  4 

TT  W 


11  C  J  4 
+  — r  it  w„ 


pi  1,4tt  2  pi  2,4tt 


(i> 


l 


+  wi.66t(wi.e  +  u2.e’ 


,} 


sin  tt  4  cos  t  d4d  t  =  0 


(95) 


3A 


Carrying  out  the  integration,  and  substituting  the  values  — —  from 

“l 

equation  (88),  and  c  q  and  e  q  from  equations  (43)  and  (44),  we  obtain 
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Thus  in  the  second  approximation,  the  response  equation  takes  the 
following  form: 

°.2  ]+  *11  2  *2. 

T  =  1  +  Ti  7  +Tl  — 7 
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(97) 


In  equation  (97),  terms  with  the  factor  —  are  due  to  rotatory  inertia 
(J  =  pi,  —  =  1,  the  factor  — —  is  retained  solely  for  identification 

Ol  pi 


purposes);  while  with  7^-7  are  due  to  shear  deformation.  All  these 

K  O  |  j  1 

terms  are  with  the  slenderness  ratio! — ~l  as  their  coefficients. 

a  a 

Since  beams  associated  with  large  deflection  are  usually  of  the 
slender  type,  for  all  practical  purposes,  the  shear  deformation  and 
rotatory  inertia  terms  can  be  neglected,  and  thus  equation  (97)  can  be 
reduced  to 
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Equation  (98)  indicates  that  the  nonlinear  effect  on  frequency  due  to 

amplitude  is  of  the  soft  spring  type  (i.e.,  the  response  curve 
bends  towards  the  left).  As  is  seen  in  figure  16,  it  also  demonstrates 
the  jump  phenomena.  By  increasing  the  frequency  it  jumps  up  from 
point  f  (on  left  branch)  to  point  d  (on  right  branch),  and  by  decreasing 
the  frequency  it  jumps  down  from  point  b  to  c.  The  portion  b  e  f  of 
the  left  branch  of  the  response  curve  is  unstable.  (Appendix  B) 

The  solution  for  w  and  u  up  to  the  second  approximation  can  be 
written  as 


A  L- 

4  w  q 

w  =  a  q  sin-n-f;  cos  fi  t  +  -  cos  fi 

l  1  rr 


co 

,  I 


sin  m 


jlJl 


m  =  3,5,. . .  m(m 


1) 


u  = 


1 


8  ‘"V’2  (6+  2„ 


sin  2tt(;)  (  1  +  cos  2fijt) 


(99) 

(100) 


For  a  given  frequency  fi^,  the  amplitude  (qa^)  is  determined  by 

equation  (98),  and  the  solution  of  u  and  w  can  be  obtained  by  equation 

(99)  and  (100).  Equation  (100)  shows  that  period  of  u  I— — M  is  only  half 

[2  fij 

of  that  of  the  external  excitation  or  of  w.  It  also  indicates  that  when 
(1  4  cos  2fijt)  -  0,  u  —  0,  and  thus  u  —  0  when  w  is  zero*  The  maximum 

of fu|  occurs  at  £  =  1  an 


1 


(TT  a  ql 

!  *  1 

1  1 

u  = 

1  1  max 

2  / 

.  Examining 

equation  (99),  it  is  seen  that  the  second  term,  i.e.,  the  term  with  the 
summation  sign,  is  very  small  and  converges  rapidly.  Therefore, 
it  can  be  neglected. 
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SECTION  4 


MODIFICATION  OF  THE  PROBLEM 
DUE  TO  MATERIAL  DAMPING 

4.1  STRESS-STRAIN  LAW  WITH  HYSTERESIS 

Under  alternating  loads,  the  deformation  of  materials  do  not 
follow  Hooke's  law.  Part  of  the  deformation  energy  is  converted  into 
the  internal  friction  loss.  Due  to  this  process  of  dissipation  of 
energy,  the  stress -strain  curve  forms  a  loop,  which  is  commonly 
known  as  the  hysteresis  loop.  The  area  of  the  loop  determines  the 
amount  of  energy  dissipated  per  unit  volume  of  the  material  during 
one  cycle  of  oscillation.  This  property  of  energy  dissipation  is  called 
the  internal  damping,  or  material  damping,  and  is  a  characteristic  of 
the  material.  Material  damping  serves  as  an  effective  agent  to 
reduce  the  amplitude  of  vibration. 

One  relationship  that  is  commonly  used  for  material  damping 
is  in  the  form  of  viscous  damping,  i.e.,  the  damping  is  proportional 
to  the  rate  of  deformation  during  oscillation.  This  is  the  same  as 
saying  that  damping  is  proportional  to  the  frequency  of  the  oscillation. 
According  to  this  assumption,  the  equations  of  oscillations  are  linear, 
and  can  be  solved  often  by  rigorous  mathematical  methods.  However, 
experimental  work  shows  that  damping  in  general  is  only  dependent  on 
the  amplitude  of  the  oscillation  and  is  independent  of  the  frequency. 
The  use  of  viscous  damping  for  the  vibration  analysis,  therefore,  is 
not  justified. 
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Based  on  the  fact  that  the  material  damping  is  dependent  on 
the  amplitude  and  is  independent  of  frequency,  Davidenkov  proposed 
an  empirical  stress -strain  law.  Davidenkov's  law  was  adopted  by 
Pisarenko  for  cantilever  beams  (both  Euler -Bernoulli  and 
Timoshenko  beams),  and  the  analysis  agrees  quite  well  with  the 
experimental  results. 

Davidenkov’s  stress -strain  law  consists  of  two  expressions. 

One  expression  is  for  the  ascending  branch  of  the  hysteresis  loop 
(loading);  the  other  for  the  descending  branch  of  the  hysteresis  loop 
(unloading).  In  functional  form,  it  is  expressed  by 

7  ‘ E  {e ' "  [ie„ + e  .n  ■ 2" 1  v]} 

(101) 

7  =  4+"K-£,n-z"'1£°1} 

or  taking  the  derivatives,  we  obtain 

>"'1] 

(102) 

i  ■  Eti  -  v<£„  -t)"'1] 

where  arrows  towards  the  right  indicate  the  ascending  branch,  i.e», 
loading  or  d£/dt  >  0  (t:time),  while  arrows  towards  the  left  indicate 
the  descending  branch,  i.e*,  unloading  or  d£/dt  <  0,  and 

E  is  the  Young’s  modulus 

£  is  the  amplitude  of  the  strain  (absolute  value) 

£  is  the  strain  at  any  instant 
cr  is  the  corresponding  stress 
v,  n  are  parameters  for  the  hysteresis  loop. 
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A  sketch  of  the  hysteresis  loop  of  equation  (101)  is  shown  in  figure  5. 


(T 


Fig.  5  -  Hysteresis  Loop  --  Davidenkov's  Stress-Strain  Law 
The  following  are  the  properties  of  the  Davidenkov's  stress- 


strain  law: 


(a)  The  damping  energy,  i.e.,  the  area  of  the  loop,  is 

independent  of  the  frequency  but  is  dependent  on 
the  strain  amplitude  £  . 

(b)  For  the  same  amplitude  of  strain,  different 
materials  have  different  loop  areas.  The  loop 
area  then  depends  on  the  two  parameters  v  and  n. 
These  two  parameters  also  determine  the  shape 
of  the  hysteresis  loop.  However,  as  shown  by 
Pisarenko  [30},  the  shape  of  the  hysteresis  loop 
has  little  effect  on  the  damping  properties  and 
therefore  if  we  take  n  to  be  the  number  2 

or  3  and  vary  v  accordingly,  the  hysteresis  loop 
thus  obtained  usually  gives  satisfactory  results 
for  damped  vibration  problems. 

(c)  The  following  conditions  of  symmetry  are  satisfied: 


df  e 


d£ 


e  =  +e 

v  o 


d£  £ 


o 
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(d)  The  starting  points  of  loading  and  unloading  are 
always  with  a  slope  equal  to  E,  the  Young* s 
modulus.  Thus 


dcr 

d? 

dcr 

d£ 


£  =  -£ 


£  =  >£ 


=  E 


=  E 


Davidenkov's  law,  as  expressed  by  equation  (101)  is  only  for 
the  case  when  the  strain  is  symmetric  about  the  origin.  For  the 
case  in  which  the  hysteresis  loop  is  not  symmetrical  about  the  origin, 
i.e.,  the  extreme  values  are  not  equal,  £  £ 

111  I  Cs 


\  * 


Fig.  6  -  Unsymmetrical  Hysteresis  Loop 

the  same  functional  relationship  of  Davidenkov's  law  and  make  the 
assumption  that  the  loop  is  symmetrical  about  the  point((£j  +  ^ ) /2 , 

(<Tj  +  The  stress -strain  law  then  takes  the  following  modified 

form: 

7‘  E<£-  1  [<£-£/-  ?<£2-£,»1> 


■  E<€+  H(Vt,n  -  £  **-£i»n]> 


(103) 
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This  modified  form  of  equation  (103)  also  possesses  the  properties 
(a)  -  (d)  of  Davidenkov's  law,  but  now  the  amplitude  of  the  strain  is 
defined  as  the  average  of  the  extreme  values,  i.e,}  *  ^j/2* 

For  most  metals,  the  parameter  n  can  be  taken  as  2  ^3oj.If  n  =  2, 
equations  (103)  become 

E{e-  :r[e2-ze1t+  \  i^2 + 

?  =  E{e+  He2-2£2e-  ?  (io4) 

In  the  following  development,  the  Davidenkov's  stress -strain 
law  in  its  modified  form  expressed  by  equations  (103)  or  (104)  will 
be  used* 

4.2  NATURE  OF  STRAIN  CYCLES  IN  BEAMS  AND  STRESS 

IN  TERMS  OF  BEAM  DISPLACEMENTS 

Before  deriving  the  equations  for  the  motion  of  beams  with 
hysteretic  material  properties,  it  is  necessary  first  to  get  a  physical 
understanding  of  the  strain  cycle  occurring  at  any  point  in  the  beam. 

For  the  sake  of  simplicity,  the  excitation  frequency  will  be 
assumed  to  be  near  the  lowest  linear  flexural  mode  frequency,  and 
just  as  in  the  undamped  case,  only  two  cases,  one  with  both  ends 
hinged  and  the  other  with  simply  supported  ends,  will  be  considered, 

4.2.1  Both  Ends  Hinged  Case 

Since  the  hysteretic  damping  effects  are  very  small  in 
comparison  with  the  elastic  effects  (for  metallic  materials)  it  will  be 
assumed  that  the  basic  strain  cycle  is  similar  to  the  one  for  an  elastic 
beam. 
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Section  3.3  gives  the  results  in  the  first  approximation 


for  an  elastic  beam  with  hinged  immovable  ends  as  follows: 


w  =  l  w  %fWjTi  =  I  a  q  sin  ir£  cos  t 


1  2  12  2  2 
u,  =  °-  u,  t  +  7  w,  t  =  T  a  tr  cos  r 

1  Z,i  Z  1  ,£  4  1  (105) 

The  axial  strain  of  the  middle  plane  is  expressed  by  equation  (47): 


£=  ui  ,4  ^  +  (uz,4+  r  wi142)"2+  '  '  ' 

The  axial  strain  at  any  point  is  obtained  from  equation  (5): 

£=  e  -  ze 

.X 

Substituting  equation  (105)  into  (47)  and  (5),  we  have 

f-  1  i  .2  z  2  . 

t  -  —  (a  t|  ir  cos  t)  +  —  a  t)n  sm  ir^  cos  t 
x  1  JL  1 


(47) 


(5) 


(106) 


where  the  first  term  is  due  to  axial  tension,  while  the  second  term 
is  due  to  bending. 

Differentiating  £  with  respect  to  t, 


dc  £  1  z  . 

—  =  *  a  n  tt  smr{r-  a^  r\  cos  t  +  —  sin  tr  £) 

dr  l  1  XZ  1  i  *  (10?) 

Positive  values  of  d€/dT  corresponds  to  the  loading  branch  of  the 
hysteresis  loop  (figure  b);  while  negative  values  of  d£/dT  corresponds 
to  the  unloading  branch.  The  number  of  strain  cycles  that  occur  at 
any  point  during  one  cycle  of  beam  motion  is  determined  by  the  number 
of  changes  in  sign  that  occurs  in  the  strain  rate  during  that  interval 
of  time  at  that  point.  The  hysteretic  stress -strain  law  for  each  strain 
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cycle  is  also  dependent  on  the  maximum  and  minimum  value  of  £ 
of  that  cycle  (£j,  £z  of  equation  (104)).  For  these  two  reasons,  the 
stationary  values  of  equation  (106)  will  now  be  examined.  The 
stationary  values  of  £  is  determined  by  the  zeros  of  equation  (107), 
i.e.,  at  t  ~  Of  it.  2ir ,  «... 

and  when 


1  z  .  _ 

—  a  ncos  t  +  —  sin  ir£  =  0 

2  1  f 


(108) 


i.e when 


t  ;  r  =  cos 
o  o 


1  /  -2  z  sin  m 

V*  I 


(109) 


for 


TT  ^  3tT 

z  >  0,  —  <  t  <  it;  and  tt  <  t  <  — 

2  o  o  2 


Since  j  cos  t|  <  1,  no  solution  of  equation  (109)  is  possible  if 


|  2z  sin  ir£  |  >  a^  I 


There  are  two  types  of  cyclic  variations  of  strain  and  they  depend 
on  the  solution  of  equation  (109).  These  two  types  are  described  below. 

(a)  Case  when  |  2  z  sin  J  <  a^ql 


For  this  case,  for  the  interval  0  <  t  <  2ir,  there  are 
two  roots  of  equation  (109).  These  roots  shall  be  called  tq  and  tq. 
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For  z  >  0,  we  have:  —  <  t  <  tt,  and  ir  <  t  <  — —  ,  Based  on  the 

c—o—  —o—2 

knowledge  of  all  the  zeros  of  equation  (107)  and  the  signs  of  df  /Ht,  it 
is  possible  to  determine  the  time  dependence  of  the  strain  from 


equation  (106),  and  the  maximum  and  minimum  values  of  strain.  Thus, 

dC  -  dC 

at  t  =  0,  — —  =  0,  and  since  for  t  =  0  ,  - —  =  +,  i.e.,  c  is  increasing,  and 

+  df  ,  r  dT 

t  ~  0  ,  —  =  i.e.,  6  is  decreasing,  t  must  be  a  maximum  at  t  =  0. 

dT  df  dr 

On  the  other  hand,  at  t  =  t  ,  again  — —  =  0,  but  for  0  <  t  <  t  ,  =  - 

o  dT  df  °  dT 

and  hence  £  is  decreasing,  and  for  t  <  t  <  tt,  — —  =  4-  and  hence  £  is 

o  dT 

increasing,  £  must  be  a  minimum  at  t  =  r  *  By  the  similar  reasoning 

£  also  must  be  a  maximum  at  r  =  tt;  and  a  minimum  at  r  =  r  .  In 

o 

summary,  for  z  >  0,  we  have 


T 

0 

T  TT 

O 

t  2  ir 

o 

d£ 

0 

0  +  0 

0+0 

dr 

(decreases)  (increases)  (decreases!  lmrrPaeo=\ 

state  of  £ 

max. 

min.  max. 

min.  max, 

For  a  given  set  of  values  of  z  and  £  f  ths  strain  cyclic  variation  is  shown 
in  figure  7(a).  It  is  seen  that  there  are  two  cyclic  variations  of  £ 
during  one  cycle  of  excitation  (U  <t  <  2  rr). 


(b)  Case  when  J  2  z  sin  tt£  |  >  <=  i 

In  this  case  there  are  no  solutions  for  equation  (109). 

It  is  observed  from  equation  (107)  that  the  following  is  true.  For 
z  >  0 


T 

0  n  2  t r 

d£ 

0  , ,  -  .  o  +  0 

(decreases)  (increases) 

dT 

state  of  £ 

max.  min.  max. 

Therefore,  there  will  be  only  one  cycle  of  variation  of  £  during  the 
interval  of  0  <  r  <  2  tt.  The  time  variation  of  £  for  this  case  is 
shown  in  figure  7(b), 
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At  a  given  instant  (i.e.,  t  =  a  constant  )  equation  (108) 
represents  a  curve  in  the  z  -  £  plane.  Sketches  of  the  curve  are 
shown  in  figures  8(a)  and  (b).  A  series  of  curves  at  different  instants 
of  time  for  the  case  when  a ^ !>h  are  shown  in  figure  8(a).  Starting 
from  outside  (at  t  =  0),  the  curve  moves  toward  the  center  of  the  beam, 
and  coincides  with  the  center  line  of  the  beam  (£  axis)  at  t  =  tt/2. 

From  t  =  tt/2  to  t  =  ir  the  curve  moves  outward  in  the  other  direction. 
From  t  =  tt  to  t  =  2ir,  the  motion  of  the  curve  is  reversed.  Figure  8(b) 
shows  the  curves  for  the  case  a^qlCh.  The  curve  again  moves  with 
time  inward  at  the  side  for  z  <  0  during  the  interval  0  <  t  <  ir/2, 
coincides  with  the  £  axis  at  t  =  W2,  moves  outward  at  the  side  for 
z  >  0  during  the  interval  ir/2  <  r  <  it,  and  then  reverses  the  motion 
during  the  interval  it  <  t  <  2ir.  The  difference  between  the  cases 
represented  by  figure  8(b)  and  figure  8(a)  is  that  in  figure  8(b),  the 
curve  starts  inside  some  part  of  the  beam  at  either  t  =  0  or  t  =  tt, 
and  there  is  a  portion  of  the  beam  (two  shaded  areas)  that  is  always 
outside  the  curve;  while  in  figure  8(a),  at  either  t  =  0,  or  t  -  it,  the 
curve  is  completely  outside  of  the  beam.  This  difference  is  signifi¬ 
cant  in  that  the  portion  of  the  beam  which  is  always  outside  the 
curve  (the  two  shaded  areas)  is  the  area  for  which  |  2z  sin  ir£|  >  a  ql. 

and  thus,  for  which  no  solutions  (t  .  t  )  of  equation  (109)  exist. 

o  o 

Consequently  the  strain  cycle  for  this  portion  of  the  beam  can  be 
represented  by  figure  7(b).  For  the  other  portion  of  the  beam  (area 
which  is  not  shaded  in  figure  8(b),  and  the  entire  beam  of  the  case 
shown  in  figure  8(a),  the  strain  has  two  cycles  which  is  represented 
by  figure  7(a). 
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Figure  7  -  Periodic  Variation  of  Strains  for  Immovable  Hinged  Ends  Case 
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(b)  Case  for  t]  jf  <  h 


Figure  8  -  Sketches  of  the  Curve 
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(a)  Two  loop  case  —  |  2  z  sin  tt  |  |  <  a^  rii,  corresponding 
to  Figure  8(a)  and  unshaded  area  of  Figure  8(b). 


-  0 ;  2  tt 


(b)  Single  loop  case  --  |2zsin7r£|>  a  corresponding 

to  shaded  area  of  Figure  8(b)* 

Figure  9  -  Hysteresis  Loop  of  Beam  Stresses 
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Based  on  the  observations  on  figure  7  and  8,  and  the 
above  discussions  on  the  two  types  of  cyclic  variation  of  strain, 
stress -strain  curves  of  the  hysteretic  type  can  be  formed  as  shown 
in  figure  9(a)  and  (b). 

The  functional  form  of  the  stress -strain  law  can  be 
formulated  by  using  equation  (104).  For  the  two  cycle  case,  for 
z  >  o  and  the  part  of  the  cycle  otq,  figure  9(a)  as  compared  with 
figure  6,  we  have, 

£  =  £  (tq)*  (min.)  and  £  =  £  (0),  (max.) 

while  for  z  >  o,  and  part  of  the  cycle  irf  : 


£,  =€(t  ),  (min.),  and  £  =£(tt),  (max.) 

1  o  c 

Note  that  £ (t  )  =  £  (t^).  From  equation  (5),  we  have 

£  (0)  =  £ (0)  -z0,  (0) 

x 

£  (it)  =  £  U)  -  Z  0,  (it) 

X 

from  (84),  (105)  and  (47),  it  can  be  shown  that 

£(0)  =  £(ir),  and  0,  (tt)  =  -  0,  (0). 

x  x 

For  brevity,  let  us  denote  £{0)  by  £  ,  and  0,  (0)  by  0,  ,  then  we  have 

m  x  xrn 


C(o)  =  £ 


z  0, 

xm 


£(tt) 


+  z  0, 

xm 


(110) 


To  obtain  the  minimum  value  £  (t  ),  substituting  equation 


(109)  into  equation  (106),  thus  we  have 


£(t  )  = 
o 


IT  sin  TT  £ 


(111) 
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This  indicates  that  £  (t  )  is  proportional  to  (z/f)2  ,  and  if  (z It)  is  small, 
€(tq)  can  be  neglected. 


Using  equations  (104),  (110)t  and  neglecting  the  £  (to)  =  £^ 
terms,  we  have,  in  terms  of  £  and  0,^,  the  following  stress-strain  law 
with  two  loops 


cr 

OT 


o 


c r 

ITT 


O 


O' 


OT 

O 


tr 

ttt 

o 


e  ( (c- z0,  )  +  ^  -  (€  2  -j  e  2)  +  zuee. 

x  c  I  £.  m  x 

-  _  2212 

+  G  0.  )  -  *  (e.  ~~  e.  )  l 

m  xm  x  2  xm  I  j 

e(  (£- *e.  )  +  ff(£2  +  7C  2-2?e  ) 

I  x  2  I  2  m  m 

-  z(2€e,  +  £  e.  -  2£  0.  +  €e,  ) 

x  -  m  xm  m  x  xm 

2  2  12- 
+  z  { 0,  +  —  0,  +2  0,  0, 

x  2  xm  x  xm 


(H2) 


(113) 


Now  we  turn  our  attention  to  the  single  loop  case  (figure  9(b)  ).  As 
figure  9(b)  is  compared  with  figure  6,  we  have 


For  z  >  0:  =  £(tt),  min. 


G2  =  €(0).  max. 


and  for  z  <  0:  =  G  (0),  min. 

€2  =  £(*■),  max. 


However,  due  to  symmetry,  we  have 

G(z)  =P(-z) 

'  at  t  =  0  ^  at  r  =  Tr 
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Thus,  for  the  part  of  the  cycle  ott,  we  have 


e  -  e  -  H  e- 

2  I  I  xm 

for  all  z 

€  =  e  +  Izl  e,  (H4) 

1  m  1  i  xm 

Substituting  equation  (114)  into  equation  (104),  the  stre  ss  -  strain  law 
of  the  single  loop  type  is  obtained: 


0 


a  =  E{  (£-  z0  )  +  -^|  +(€  -  e  )2±2z(e-€  )0. 

i  v  vi  m  '-m  x 


+  2|z|  0,  (£-£  )  +  (0,2  -  0,  2  )  z2  -  29,  0,  z 

xm  m  x  xm  x  xm 


H]} 


(115) 


4.2.2  Simply  Supported  Case 

The  first  approximate  solution  for  the  simply  supported 
case,  as  given  in  Section  3.4,  shows  that 

w  =  Iw  Z  iw ^  t|  =  v\  sin  tt|  cos  t 


"i  =  °-  u2,e  +  l”i,e2  =  0 


(116) 


Substituting  equation  (116)  into  equations  (47);  (5),  we  have, 

r  z  2  . 

fc  =  —  it  a^  q  sin  ir£  cost 


d£  _  z  2  . 

—  =  _  —  ir  a  n  sin  wf  sin  r 

dT  i  1 


(117) 


Equation  (117)  indicates  that  the  strain  varies  linearly 
across  the  depth  of  the  beam  and  the  neutral  axis  is  at  the  center.  It 
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also  forms  a  single  loop  from  0  to  tt  and  back  from  tt  to  2w,  and  the 
loop  is  similar  to  the  one  shown  in  figure  9(b).  The  functional  form 
also  must  be  the  same.  In  other  words,  for  the  simply  supported  case 
with  one  end  axially  free  to  move,  the  stress -strain  law  is  ex¬ 
pressed  by  equation  (115). 

Note  that  the  main  parts  (elastic  part)  of  the  stress- 
strain  law  of  equation  (112),  (113),  and  (115)  are  all  the  6ame  a s 
expressed  by  the  term  E  (G-  z0>x)-  The  onlY  deference  is  the  second 
order  effect  of  hysteresis  which  is  represented  by  the  terms  with  the 
coefficient  v. 


4.3 


AXIAL  FORCE  AND  BENDING  MOMENT  TAKING  MATERIAL 
DAMPING  INTO  ACCOUNT 


The  axial  tension  and  bending  moment  are  defined  respectively: 


N  = 

M  = 


u  d  A 


(r  z  d  A 


(118) 


The  integration  is  over  the  entire  cross  section  of  the  beam.  Again 
the  immovable  hinged  ends  case  and  simply  supported  case  shall  be 
considered  separately. 


4.3.1  Both  Ends  Hinged  Case 

In  Section  4.2.1,  it  is  shown  that  the  addition  of  the  axial 
tension  causes  the  beam  to  have  at  some  points  two  cycles  of  strain 
and  some  points  one  cycle  of  strain  during  one  cycle  of  beam  motion. 
The  criterion  of  determining  which  part  of  the  beam  has  two  cycles 
and  which  part  has  one  cycle  is,  as  explained  before,  obtained  by 
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examining  the  roots  of  equation  (109);  while  the  criterion  of  deciding 
which  branch  of  the  hysteresis  loop  is  applicable  is  to  examine  the 
signs  of  the  loading  rate  (d€/di-)  of  equation  (107)*  In  both  equations, 

(109)  and  (107),  there  are  three  variables  —  t,£  and  z,  and  thus  at 
any  given  instant  r,  but  at  different  locations  of  beam  (different  £  and 
z),  stress -strain  relationships  may  be  different*  Consequently,  the 
integration  of  equation  (118)  will  not  only  give  complicated  expressions 
for  N  and  M  but  also  take  different  forms  at  different  part  (§)  of  the 
beam*  In  general,  therefore,  the  analytical  expressions  for  the  damped 
motion  of  the  two  hinged  ends  case  are  very  complicated*  For  the  sake 
of  simplicity,  only  two  limiting  cases  will  be  considered.  One  limiting 
case  will  be  for  beam  amplitude  of  vibration  large  in  comparison  with 
beam  depth;  while  the  other  limiting  case  is  the  one  for  beam  amplitude 
small  in  comparison  with  beam  depth  but  still  large  enough  to  produce 
nonlinear  effects  (figure  7(b)  ).  It  may  be  expected  that  the  general 
case  (intermediate  case)  will  have  a  behavior  intermediate  between 
these  two  limiting  cases. 

4*3.11  When  Amplitude  is  Large  Compared  with  Depth 
of  the  Beam'" 

If  the  amplitude  is  sufficiently  large  so  that  in  equation  ( 109) 
2  2 

cos  T  =  — -  sin  ir£  =  0  {i.e.  /a  n>>2z  sin  *£) 

/an  1 

1 

then  t  ~  rr  =4"  (figure  7(a)  ). 

o  2  o  2 


In  the  following,  for  brevity,  this  case  will  be  referred  to  simply  as  the 
very  large  amplitude  case 
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If  the  approximations  tq  =  j  and  tq  =  ~  are  used,  the  unloading  and 


IT 

Z  O  2 

the  loading  cycles  of  the  beam  can  be  written  as 


0  <  T  <  T* 
2 


a  “  <r 


0T 


TT 

2  <  T  <  " 


tr  -  tr 


TTT 


TT  <  T  < 


3w 


cr  =  or 


TT  T 


3t r 


<  T  <  2  TT 


cr  =  cr 


Ot 


U19) 


Using  the  appropriate  formula  of  equations  (112)  and 
(113)  and  substituting  (119)  into  equation  (118),  we  obtain  the  following 
results. 


IT 


3^  \=  EAt  EA  (C2  -it2  )  -  |EI(0,2  -ie,2  ) 

N  (—  -  Zir,  '  2  2  m  2  x  2  xm 


N  (0  -  — )  1 

>=  eae  t-j  ea  (€  +  k2  -  zii  )  + 1  ei 

N  ,  2  2  m  m  2 

2  ^  2  12- 

(0,  +-0,  4-20,  0,  ) 

x  2  xm  x  xm 


(120) 


M  (-  ~  tt) 

3ir  )  =  -Eie  +^Ei(2ee,  +  G  e,  ) 
M  (^-  -  2w)  '  *  2 


m  xm 


M  (0  -  j) 

)  =  -Ei0,  --^Ei{2ee,  +£0,  -  2C  e,  +  €e,  ) 

.  .  ,  (  x  2  x  -  m  xm  m  x  xm 

M  (ir  _  — ) 


(121) 
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where 


4. 3. 1.2  When  Amplitude  Is  Small  Compared  with 

Depth  of  the  Beam" 

If  the  amplitude  is  small  so  that  we  make  the 
approximation  2z  simr£  >  t  a^q,  then  there  is  no  solution  for  equation 
(109),  and  the  stress-strain  law  takes  the  simpler  form  of  equation 
(115).  Strictly  speaking,  however,  the  relationship  of  2z  sinir|  >  la^ 
can  not  hold  for  the  whole  span  of  the  beam,  since  simr£  approaches 
zero  toward  the  end  of  the  beam  (at  £  =  0  and  1 ),  while  f  a^q  is  never 
zero,  no  matter  how  small  it  may  be.  The  approximation  is  taken  in 
the  sense  that  l  a^q  is  sufficiently  small  so  that  at  any  instant  only  a 
very  small  portion  of  the  entire  beam  behaves  with  the  hysteretic 
behavior  of  equations  (112)  and  (113)  rather  than  with  the  hysteretic 
behavior  of  equation  (115),  and  the  substitution  of  equation  (115)  for 
equations  (112)  and  (113)  for  this  small  portion  of  the  beam  does  not 
alter  the  damping  behavior  too  much.  Physically  speaking,  this  is  the 
case  when  the  shaded  areas  of  the  beam  shown  in  figure  8(b)  become 
predominant;  while  the  unshaded  area  becomes  negligibly  small. 


:,:In  the  following,  for  brevity,  this  case  will  be  referred  to  simply 
as  the  small  amplitude  case. 
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Based  on  this  approximation,  and  the  loading 
cycle  as  shown  in  figure  9(b),  the  substitution  of  equation  (115)  into 
equation  (118)  gives 


n  =  ea€+  eiv  (C - e  )  (0.  +  9.  ) 

1  m  xm  x 


M  =  -  E 1 9,  -  El  v  9,  9,  + 

x  3  x  xm 


2  2  k 

+  I,  (9.  -  0.  ) 

3  x  xm 


0 


7  ['■ 


(C'C  J 

m 


(122) 


where  upper  signs  refer  to  n  <  t  <_  2tr  ,  lower  signs  refer  to  0  <_  t  <  tt, 
and 


I 


1 


I 


3 


(123) 


4.3.2  Simply  Supported  Case 

As  stated  in  Section  4.2.2,  for  the  simply  supported  case 
the  stress- strain  law  follows  the  same  form  of  equation  (115),  which 
in  turn  gives  the  axial  tension  N  and  bending  moment  M  the  same 
expressions  as  shown  by  equation  (122). 

4.4  EQUATIONS  OF  MOTION  WITH  MATERIAL  DAMPING 

As  observed  in  the  analysis  of  Part  I,  for  the  type  of  slender 
beams  considered,  the  influence  of  the  shear  deformation  and  rotatory 
inertia  is  small.  In  the  following  analysis  the  shear  deformation  and 
rotatory  inertia  will  be  neglected  for  simplicity. 
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Since  shear  deformation  is  neglected,  we  have 


y  =  0,  and  0  =  a 


(124) 


and  equation  (7)  becomes 

w, 

tan  0  =  - - - — 

1  +  u, 

x 


sin  0 

cos  0 


w, 

x 


i  +  e 


\  +  U.x 

i  +  e 


(125) 


When  both  shear  deformation  and  rotatory  inertia  are  neglected, 
equation  (12)  is  reduced  to  Q  =  M,x/(1  +  £),  By  using  this  expression 
and  equation  (125),  Q  and  £  are  eliminated  from  equations  (10)  and  (11), 
and  thus  we  have 


,  (M, 

,mT  1  I  x  sin  2 0  .  . 

(N  cos  0),-=  l — — - 1  =  pAu, 

x  2  V  1  +  u(  /  tt 

'  x  ,x 


(126) 


(N  sin  0),  4- 

x 


M(  cos  0  . 

- - - j  +  K(x,t)  =  pA  (wt  -  £7  w  cos  £7  t)  (127) 

i  +  up  / ,  tt  i  l  ;  ' 


X  /  X 


Since  there  will  be  a  phase  difference  between  the  external 
excitation  and  the  beam  motion,  a  phase  angle  is  introduced  as 
follows 


r  =  t  +  iJjj  ,  or  =  t  - 


(128) 


and  we  let  be  expressed  by  a  perturbation  series: 


+1  =  +01  +  +11"  +  +  • 


(129) 
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Then 


cos  t  =  cos  (t  -  4^) 


=  cos  (t  -  4<01).+  n  sin  (r  -  4iQ1) 

-  n2  cos  (T  "  4>01)  +  «1>21  sin  (t  -  4^)1  +  •  •  •  U30) 


Now  values  of  N  and  M  are  substituted  in  equations  (126)  and  (127)  and 
the  resulting  equations  are  transformed  into  dimensionless  form. 
Perturbation  equations  are  obtained  by  techniques  similar  to  those  used 
in  Part  I  for  the  undamped  case.  Among  the  difference  in  the  resulting 
equations  are  those  due  to  the  introduction  of  the  phase  angle  4^  >  and 
those  due  to  the  fact  that  the  expressions  for  N  and  M  at  different 
intervals  of  time  (t)  are  different  now.  Consequently  the  equations  of 
motion  take  different  forms  during  different  time  (t)  intervals.  After 
some  computation,  we  have 


ui.££  =  ° 


4 

0  w  +  w .  „  t  „  =  0 

pi  i  >  tt  i.eeee 


i 

u„  _  +  —  w  . 


2,e  2  i.i  j'e 

4 


,,=  o 


cPl  U2,tt  -  ,U3.t  +  Wl.eW2.etCWl.iWl.«4,'t 


(131) 

(132) 

(133) 

(134) 


where  the  notation  is  the  same  as  in  Part  I. 

Equations  (131)  (132),  (  133)  and  (134)  are  valid  for  all  intervals 
of  time  (0  <  t  <  2ir),  and  for  all  boundary  conditions.  In  fact,  they  are 
identical  to  the  corresponding  equations  for  undamped  motion, 
(equations  (51)  (52),  (  53)  and  (55)).  The  second  order  and  third  order 
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equations  for  w  {i -  e..  w^,  and  w^,  however,  are  all  different  and  are 
listed  below. 

(a)  For  both  ends  hinged  case  with  very  large  amplitude 


+  w2.6eee +  pi“ 

Z  1 1  T  T 

“i 


7<Iwl.6twl.6u2.6,.6 


-Pl4»o  cos  (r  -  +01,  +  ^  [  l*wltK  ♦  w1>Mm) 

,  ,  1  2  1  2  1 

(u?  A  4  7  Wi  A  -  u->  t  ~  T  w  )  =  0 

2,5  2  1,5  2,5m  2  1,5m  I  ,£| 


(135) 


where  the  upper  sign  la  for  0  <  T  <  ~  ,  and  the  lower  sign  is  for 


3ir 

Tr  <  t  <  —  and 


pi4w2.tt+"'2.5U5+p,4  An 


1.1  3 


uj 


1 


—  w  -  —  (—  w  .  +  w  u  y 

2  1 , tt  C  *2  1,5  1,5  2,5',5 


p14Woco.(r-  V.  U  [2wu|(uu4,,;, 

2  W1.S  *m]  ,46  =  ° 


+  W1  tt  (u,  t  + 

—  1,55m  2,5 


(136) 


where  the  upper  sign  is  for  j  <  r  <  it,  and  the  lower  sign  is 
,  3ir 

for  —  <  t  <  2ir. 


(b)  For  both  ends  hinged  case  with  small  amplitude  and 
also  for  a  simply  supported  case 
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H- 1 


P,  w  +  w 
1  2,tt 


4  An 


1  ,1 


2, iiii  +  P1  w  2  W1,TT  '  c  (2  WI,|  +  Wl,|  U2,^'^ 


2  2 
“  w„  _  ). 


4 

-  WQ  P1  COS  (t  -  4iqi)  =  0 


(137) 


.  ,  A  A 

«  4  ,  „  4  11  4  21 

P,  w  +  w  +  6  - -  w  +  S  - -  w 

1  3,tt  3, an  *  2  2,tt  P1  2  I.tt 

W1  “l 

■  c  <wi,|2'“2,{  +  u3,i  wi,e>'{  t<wi,{wi,e4>’4 
2K  <W1.«  W2,H  -  wi.eemw2,eem)'ii +  2K  «**.«  wi 


,  .  11  „  4 

wue  w2,eem’’a  2  pi  W0C0,<T  -  V 

"1 


P1  "o  +n  sin(T  -  +01*  =0 


(138) 


I  V 

u  3  3  h  v 

where  K  =  =  YbT  '  and  the  uPper  si8ns  are  for  0  <  t  <  tt,  and  the 

lower  signs  for  tt  <  t  <  2tt. 
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SECTION  5 


SOLUTION  OF  THE  PROBLEM 


5. 1  BOTH  ENDS  HINGED  CASE 


5.1.1  Large  Amplitude  Case 

Up  to  the  first  three  equations  (131,  132  and  133),  the 
solution  is  the  same  as  that  of  the  undamped  motion,  and  thus  we  have 


“l  *  ° 


w  =  a  sin  tt£  cos  t 
1  1 


pi  =  " 

1  2  12  2  2 

u->t+Twie  =  T  a,  w  cos  T 
Z.%  2  1,|  4  1 


-1  2  .  ,  t  2 
u  ~  —  a  rr  sin  Znt,  cos  t 
2  8  1 


(61) 

(65) 


(66) 

(77) 


To  obtain  the  periodic  solution  of  w  from  equations  (135) 
and  (136),  first  we  must  eliminate  the  secular  terms  in  equations  (135) 
and  (136).  With  phase  shift  added,  these  terms  consist  of  both  cos  t 
sin  tt£  and  sin  r  sin  tt£.  These  terms  can  be  eliminated  by  first 
multiplying  equations  (135)  and  (136)  with  cos  t  sin  and  secondly  by 
sin  t  sin  t r£,  and  then  integrating  the  resulting  two  equations  throughout 
the  whole  beam  length  (0  <  %  <  1),  and  over  the  period  of  one  cycle 
(0  <  t  <  2ir),  which  gives: 
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,1  fz„ 

I  sin  ir£d£  / 
o  o 


n  4  ±  ^  .  4  AU 

P1  W2.tt  +  w2  .i«i  +  P, 

w 

1 


c  ( 5  Wl,?3  +  wl.{  u2,4  )  -  Pl“  »o  C°3  |T  •  V 

r  b 


— "  Ztt 


COS  T  dT 


/ 


3tt 

Z 


2w  .  I  u 

i.u  '  2,e 


i  r 

1WU2)  ]  uCOSTdT}‘ ?  /  sin"|d6{  / 


TT 

-2 


3ir_ 


{  j'[wi.«  h.e*2wi.iZ|  ] 

TT  ^  m  J 


COS  T  dT 


■} 


it  3tt 

Y 


T*  /  sin  'Sd5{  [  /  "  / 

O  ^  L  o  TT 


LWl,Htn  \U2,t 


*  2  WI.|2)  ]  aC°aTdT}tf'  j  8i"  (  j  -  J  -  J 


TT 

z 


3tt 


1 

2 


2it 


[■ 


i,e$m  |u2,e +  2  wi,e 


] 

COS  T  dT  J 

*  m 

.«  J 

}  =  0 


(139) 
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and 


/  sin  ”|d*  /  [pi4  *2,„ *  w2.sss +  pi 


4  A11 


w 


CO 


2  1,  TT 


i(lwi.|3  +  ”z.«wi.e )  -pi 

*  b 


w  cos  (r 
o 


3tt 


v] si 


sin  t  dT 


/*  r  r  2  m  "2™  r 2lT  i  r 

—K[/  -/  ♦/  •  J  ][2wi.ssK 


1  2  1 

2  wi,e  1 


3ir 


TT 

-2 


sin  t  dr 


r} - f 'I sin ,,5ds{  / 


f 

w 

?/' 


W 


us  Ks  +  iwu2  lm]  uatoTdT} 


sin  tr£d£ 


lU2,e 


1  —  TT  — 

+T  /  si"  "SdS  {  J  -J  -  J 


/2tr  “i 

■ 


_3jf 

T 


l.SSm  |U2.Stl"1,S2) 


m 


sin  t  dT  V  = 


•ee 


=  0 


(140) 
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where,  for  the  sake  of  brevity,  integrations  with  same  integrand 
but  different  integration  limits  are  grouped  together.  Using  the 
relationships  that 


/  sin  J  (  |5[ 

-*r\  n 


W2.tt  +  W2  .tUt 


V  COS  T 

)  dr  =  0, 

'  sin  t 


and  substituting  all  the  known  quantities  into  equations  (139)  and 
(140),  and  carrying  out  the  integration,  we  obtain 


5  a  ,  3  5 

a  tt  A  ,  3a  tt 

1  11  1 


+ 


2o>. 


32c 


-  2tt  w  cos  ib  -  vn  a  tt 
o  o  1  1 


3  6  3tt  .  8 


48 


=  0  (141) 


3  6 

4  .  ai  " 

2ir  w  sin  ib  -  vq -  =  0 

o  o  1  24 


(142) 


Solving  these  two  equations  for  A  and  i|j  ,  we  obtain 

11  o  1 


ij;  =  sin 

o  1 


3  Z\ 
vqa  w 

48  w 


(143) 


and 

A 


11  9*  2  (3tt-8)tt  2-  4 

—  =  a  q-- — — - vqa  +  - 

2  1  24  1  tt  a 


w 


1 


4h 


1 


w 


2  3.2 

vq  tt  a  1 

48  I 


l 

12 


(144) 


The  response  equation  is 
2 

“l  ,  ,  9JfZ  2  2  (3tt-8)tt  2  2  7  4 

- —  =  1  + a  q - —L — vq  a  +  - 

1  24  1  tt  a  q 


to 


1 


4h 


1 


2  3  3  21 

i  vtt  a .  q  > 

2  4  | 

Wo  ’I 

1  48  | 

(145) 


Note  that  if  v  =  0,  equation  (14  5)  is  reduced  to  the  response  equation 
for  the  undamped  case,  i.e.,  equation  (76),  and  that  because  of  the 
damping,  the  response  curves  are  closed  at  the  top  now.  To  solve 
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for  we  observe  that  equations  (135)  and  (136)  both  can  be  put  into 
the  following  form: 


4 

7T  W, 

2,  TT 


+  W 


F(t,£) 


(146) 


By  expressing  F(t,£)  into  a  double  Fouries's  series  for  the  intervals 
(0  <  I  5  U  and  0  <  t  <  2tt,  these  two  equations  can  be  combined  into 
one  equation.  Thus 


4 

TT  W„ 

2,  TT 


+  W 


2.H££ 


4tt  w  cos  (t  -  ^  ,) 
O  ol 


Z 


sin  rmr£ 
m 


1  3  4,  vr|  3  5  . 

T6eE!  ^  Sin  ^  C°S  3t  +  ~  ai  sm 


oo 


z 


m  =  3.  5, 


m  +  1  2 

.  ..  2  m  -  3 

(-  1)  - x -  cos  mT  + 


oo 


m  (m  -  4) 


z^ 


2  4  2 

l)(m  -  7m  +  18)  sin  mT 


m  =  5.  7, 


2  2  2 
m  (m  -  4)(m  -  l)(m  -9) 


22  •  o 

t  -  sm  3t 


(147) 


0  <  t  <  2tt,  0  <  £  <  1 


where  the  secular  terms  are  being  removed  by  the  use  of  equations 
(141)  and  (142).  Solution  of  equation  (147)  is 
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oo 


w  =  —  W  COB  (  T  -  + 

2  11  O  O 


L 

m  =  3,  5, 


Bin  3  sin  ir£  cos  3t 

m,m4-l)+ai  ^ 


r  00 

+  a  ^  it  sin  ir£  I  2  ^  ' 

L  ™  a  3 


m  =  3,  5, 


{m  -  3)  cos  mr 
m  (m^  -  4)  (1  -  m^) 


CD 


I 


m  ■*  3,  5 f  < 


2  4  2 

(m  4  1)  (m  -  7m  4  18)  sin  m  t  22 

T~2  m  „2  ,  2  “  '  120 

m(m  -  4)  (m  -  1)  (m  -  9) 


sin  3t 


(148) 


If  all  the  terms  higher  than  m  =  3  are  neglected,  equation  (148)  becomes 


1  ,  .  ,  .  „  ,  3  cos  3t  . 

=  7—  w  cos  (t  -  4*  , )  sin  3*£  +  a.  To 6  —  •  am 
60*  o  ol  l  128c 


w. 


WTl  3  .  ,  11.-. 

—  a,  w  sin  ir£  (cos  3r  +  71m  3t) 
40  1  b 


(149) 


Equations  (148)  or  (149)  satisfies  the  boundary  conditions  of  (58). 

The  solutions  of  w  and  u",  up  to  the  second  approximation 
can  be  written  as 


w 


1  2 

=  a,n  sin  cos  (t  -  +  777"  w  H  sin  cos  (  t  -  ^ ) 

1  dOtt  o 


33  2 

ai  ^Ai 
1281 


VTT  3  3  11 

sin  it£  cos  3t  -  —  a,  sin  (cos  3t  +  -7-  sin  3t) 


40  1 


u  =  -  -72"  (a  ti)^  sin  2*£  (1  +  cos  2t) 

1  Q  1 


(150) 

(151) 


66 


Compared  with  the  undamped  solution,  i.e.,  (79)  and 
(80),  it  is  seen  that  the  expression  of  u  remains  the  same  while  the 
expression  of  w  is  changed.  The  change  is  in  the  form  of  the  phase 
angle  4*  and  some  additional  terms  with  the  damping  parameter  v . 


Wl- 


5.1.2 


Small  Amplitude  Case 

Again  we  can  use  the  undamped  solutions  for  u  *  u  ,  and 

X  £0 


u.  =  0 


w 1  =  a  j  sin  m£  cos  t 


p,  -  » 

1  2  12  2  2 
u2.i+iwi.s  ";ai  ’  c°*  T 


-1  2  .  .  2 

u  =  —  a  tt  sin  2*5  cos  t 
2  8  1 


(61) 

(65) 


(66) 

(77) 


To  eliminate  secular  terms  in  equation  (137),  and  at  the 

2 

same  time  to  seek  the  solution  of  A  ,  ,/<•>  and  4*  ,»  the  same  technique 

11  1  ol 

used  in  Section  5.1.1  will  be  applied  again.  First  equation  (137)  is 
multiplied  by  sin  m£  cos  t  and  sin  tt|  sin  t  and  then  the  two  resulting 
equations  are  integrated  over  t  (0  -  2m)  and  4  (0  -  1).  Thus  we  have 


1  2ir 

I  sin  m£d£  |  / 

Jo  Jo 


_  i  _  4  A11 

P1  W2,tt  +  W2,£IU  +  P1  <^W1,TT 


1  ,1  3  ,  ,  , 

-  —  (—  w ,  L  +  W,  ,  u„  t)  +  2K  w,  tt  w,  t£  ) 

c  '2  1,4  1,4  2,4'  ^  1,44  1,44m'  ^ 
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2 


4 

-  W  P  cos  (t  -  ijj 

o  1 

-  cos  TdTj  3  0 

and 


)  COS  TdT 


(152) 


I 


2" 


sin  "4d£ 


p  w  +  w. 


1 3 
T  <T  w 


i  2,tt  2,4444  i  w2  i,tt  c '2  1,4 

1 


+  w,  >  u,  c)  c  +  2  K  (w  w  )  ££  -  w  p  cos  (r  -  ^  ,) 

l.e  24  >4  1.55  l»55m  >ti  °  1  oi 


r  TT  Zlt 


O  TT 


Z  2 

(w,  t£  '  w,  t£  )  tt  sin  TdT  f  = 

1.44  1,44m  ,44 


sin  t dr 


(153) 


The  evaluation  of  the  integrals  gives  the  following  two  equations 


‘l"  4H  iF5  !  8K  i  l  ,  4 

- - —  +ulal  -T1!  "  -2wo"  co<,+ol 

2 


=  0 


(154) 


32  2  5  4 

—  K  a  tt  -  2tt  w  sin  iL  .  =  0 
9  1  o  ol 


(155) 


Solving  for  A  /w  and  sin  ^  and  letting 
111  o  1 


K 


16, 


z  z 

At  12 1 


We  have 


virha^  l 

3V  ) 


(156) 
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A  .2 

H  9*  ,  ,2  h  _  —  4 

n  - *  - -  (a  Ti)  -  —  va  Tpr  + 


2  2  w  i  “r  ira  n 

4h  1  ' 


Z  4  vh  2  2  Z|  i 

*o  ’i  '  (Tiai  "  ’) 


r 


(157) 


The  response  equation  takes  the  following  form 


B. 


CD 


=  1  +  ~r  (a  t#2  -  y  va  nir  + 

,,2  1  f  1  ira 

4h 


4  f  2  4  vh  2  2  2 1  | 

-  \ w  n  -  (—  ira,  ti  ) 

Ta  t|  l  o  1  v3|  1  ]/  J 


(158) 


Again  we  see  that  if  there  is  no  damping,  i.e.,  v  -  0,  equation  (158)  is 
reduced  to  the  response  equation  of  the  undamped  case,  i,e,t  (76)*  The 
addition  of  the  damping,  however,  makes  the  response  curves  closed  at 
top  (figure  12). 

To  solve  equation  (137)  for  w  ,  first  all  the  known 
quantities  of  equation  (137)  are  expressed  in  a  double  Fourier's 
series  for  the  intervals  0  <  t<  1^  0<  <  2  tt 

Thus 

oo 

4  \  16K  m  25,  £  ,  *  3 

ttw  +  w  -  )  - — — —  a,  n  sm  rmrg  cos  t  +  4  w  ir 

2,tt  2,££g£  /  2  1  o 

™  -  *  *  m  ~4 
m  —  j  t  d  t  •  •  * 


COS  (  T  -  41  _  ) 

ol 


oo 


I 


sin  m 


m 


16c  1 


3  4  .  . 

a ,  tt  sin  cos  3r 


m  =  3,  5, 


+ 


i 


m  -  3,  5, 


16  K  sin  mT 
m(m^  -  4) 


2  5 


cos  2tt| 
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64  2  4. 

—  K  a.  j  if  stn  t 


oo 


E 


m  sin  mirf  A 

- - - >  0<  T  <  2tt,  0<  4  <  1  (159) 

m  -  4  * 


m  =  3,  5,  . 

where  the  secular  terms  are  being  removed  by  the  use  of  (154)  and 
(155). 


Solution  of  (159)  is 


W2 


QD 


m  =  3,  5, 


2  16  m  sin  mirfe  cos  t 

a  ir  K  - - - 1 - 

(m  -  1)  (m  -  4) 


4 

+  —  w 
■nr  o 


oo 


I 


sin 


m1T£ 


m  =  3,  5, 


m  (m  -  1) 


COS  (T  -  4 1  ) 

ol 


13. 

r  ,  a,  sm  if  cos  3t 
128c  1  b 


OO 


E 


m  =  3,  5, 


,  .  ,  2  sin  mT  cos  2ir£ 

. - r- - j-S - 

m  (m  -  4)  (m  -  16) 


64  3 

—  K  Sj  sin  t 


oo 


E 


m  sin  miT^ 

(m2  -  4)  (m4  -  1) 


(160) 


m  =  3,  5,  ... 

In  equation  (160)  if  all  terms  higher  than  m  =  3  are  neglected,  the 

equation  can  be  reduced  to; 

_  w 

3ir  2  o 

w  =  -  —  K  a  sin  3ir£  cos  t  +  sin  3^^  cos  (t  -  4*  j) 


70 


„  sin  wf  3  16"  2  . 

+  ^28Tai  cos  3T  -  ril  Kai  8,n3Tcos2^ 


4  2  ,  . 

—  K  sin  t  sin  3"4 


U61) 


Now  the  solution  of  w  and  uf  up  to  the  second  approximation,  can  be 
written  as 


—  ,  ,  .  .  2  sm  3ir 

w  =  a  rj  sin  cos  (x  -  tp)  +  w  r|  — r— —  cos  (r  -  <\>) 
1  o  6  Off 


A  f  3  3 


h  2  2 


1281 


a  j  t|  sin  "4  cos  3t  -  —  cos  2"4  sin  3t 


9"  3 

-  sin  3ir$  cos  T  +  ^5^  sin  3ir*  9in  T> 

—  IT  2  2 

u  =  -  —  a  r)  sin  2tt£  {1  +  cos  2t) 

16  1 


(162) 

(163) 


Compared  with  the  undamped  solution  (79)  and  (80),  it  is 
seen  that  the  expression  of  u  remains  the  same,  while  w  is  changed. 
The  change  is  in  the  form  of  the  phase  angle  ijj  and  some  additional 
terms  with  the  damping  parameter  v. 

5.2  SIMPLY  SUPPORTED  CASE 


The  solution  of  u^,  w  ,  u^,  are  the  same  as  that  for  the  undamped 
case,  i.e., 


u  =0 
1 


w  i  =  a  |  sin  cos 

1  2  „ 
u  +  —  w,  .  =0 

2,4  2  1,4 


12  2  1 

u2  =  -  -  at  tt  (4  +  sin  2"4)  ( 1  +  cos  2t) 


(84) 

(85) 

(86) 


71 


To  obtain  A  / o>  ,  ^  ,  and  w  ,  again  equation  (137) 
111  ol  2 

will  be  used.  In  fact,  all  the  results  for  A.  ./to  ,  and  w  of 

11  1  ol  2 

Section  5.1.2  can  be  directly  applied  here.  But  this  time,  we  have 


U24  +2WU2  =  0 


Thus  instead  of  equation  (157),  we  have 


11 


to. 


—  v  a  .tiir  + - 

i  1'  ira  t, 


24  vh  222 

wo  ^  -  37*1  ^ 


(164) 


The  solution  for  the  phase  angle,  however,  is  the  same  as  (156), 


Ui  ,  =  sin 

ol 


VTiha 

3w  | 
o 


(165) 


Similarly  solution  of-w?  (161)  becomes; 

M 


3K  2  1 

w?  -  -  a  ^  it  sin  3tt|  cos  t  +  w^  sin  3w|  cos  (t-  ^  , ) 


60tt  o 


ol 


(166) 


16"  2  .  4K  2 

K  a  cos  2"£  sin  3t  -  —  a  sin  3"£  sin  t 
7 ' j  1  Zb  1 


To  obtain  the  second  order  approximate  solution  of 
response  relationship  (A^j,  ‘Pjj)*  first  equation  (134)  is  solved  for 
u^.  Thus 


c  2  6  13 

U3°2al  w  1  6" 


8tt 


1  c  2  2  4 

—  sin  2ir£)  cos  2t  -  —  (2it  -  1)  a^  tr  4  cos  2t 
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9K  3  2  2  -  cost 

—  -  a  it  cos  t-  a  w  - 

50  1  1  o  40 


c.  (T ( 


16K  3  2 

+  yyg-  it  cos  t  cos  3t  (cos  ir£  + 


cos  3tt^  9 tt  3 


c  _  K  a  ,  cos  t  sin  t 
i>U  I 


64  K  3  2 


a.“  sin  )tco,  t.  f^Ka  /  cos  T  ai„T(  +  SSUIL 

1  25  1  1  2  4  1 


2925 


+  c  a ,  11 "  cos“t  (-^- +  ~-i- Sin  2ir^j 


1 


2  4" 


(167) 


In  solving  u^,  the  boundary  conditions  (81)  and  (83)  are  used. 

After  all  the  known  quantities  are  substituted  into  equation 

(138),  and  the  same  technique  as  used  in  obtaining  A^/u  2  and 

+oi  is  aPPlied  to  the  resulting  equation,  we  get  two  secular  equations 
2 

for  A  / to  and  ij;  . 

^  1  X  IX 


7  I  V 

8tt  -45  3  7  .  3  .  4tt 


-  -  .  -  ,  -  .  ,  32(997tt2-936)  V  3  5 

i92~  ai  ^  IT  (75  a1WoSin^oI  9  x  13000  T~~  *1  " 

25  ai  Wo 


a  a 
1  5  21 


C°S  t)l’  '  2  *  w  2  "  ZlT  Wo  COS  +ol  TT  +  ^V^ll  8in  +0i  =  0  <*68} 

1  1 


and 


w 


V  ,  72vI3  3  6  8wo  4 

TT  (  1 251  i  ai  "  "735*1*  C°S  ^ol  '  25  aiTr"  Sin 


o  5 


2tt  w  sin  d; 

o  o  1  (jj  2 

1 


H  4 

2 it  w  dj  cos  A  =0 
o  T  1 1  o  1 


(169) 
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Solving  (168)  and  (169)  for  A 


^  j  and  4*  .  we  obtain 


21  ~ 


CD 


8 tt “ 4 2  2  4Wo 

it  a 


96 


A  2 

1  1  +  997tt  -  936 


3hv  I 

-  8  w 

c 

81  \ 

675tt 

54hv 

2 

a  it 

125  1 

1 

1  na  cos  2  13000 

1  O  1  <a> 

1 

Z  w  2  w 


hva 


o  o 

-  -  -  COS  \h 

cos  dj  Z 5  o  1 

ol 


ol 


(170) 


+ 


1 1 


CD 


1 1 

T  tan*oi 
'1 


3vh 


8  w  i  cos  ih  , 
o  xol 


4  w 

o 

-  a  cos  iL 

135  1 


w 


81 


+  —  a, it  sin  i li  J  +  — -  „ 

50  1  Yol|  2000  w  cos  di  1 

o  o  1 


3  2  ,  vh  2 

»  (y) 


(171) 


Terms  inside  the  parentheses  of  equations  (170)  and  (171)  are  small  when 

compared  with  other  terms,  and  hence  for  simplicity  they  shall  be 

2 

neglected.  Also  the  term  A  ^  ^  is  substituted  for  the  use  of 

equation  ( 1  64).  Thus  equations  (170)  and  (171)  become 


A  Z  ? 

ZiJSn  -45)/. 

ZS  96 
o> 

1 


2  2/2wo 

i 71 1-v 


2  I 

i  A  ± 

3 vhir 

'  J. 

.2 

2  2  ' 

997tt2-  936  jh^a.)2 
+  13000  1  i  | 


(172) 
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and 


11 


vha 

TT 


(4  +_ 


,  ,  2  2 

3vhir  a^ 


Vi 


(173J 


(3w^l  )  2  -  (virha^2)  2 


Thus  the  response  equation,  up  to  terms  of  second  order,  becomes 

°1  ,  A11  A21  2 
- _  =  1  +  - _  q  +  —  *1 


w_ 


1 


u> 


1 


W. 


1 


(Sir  -45)tt  2  2  W<>T1 

=  1-  %  ain - |  l  + 


W  q  r 

— 
ira^  [ 


12  /w  2q* 
O 


na 


7  2  Z  Z  4 

(3  w  r)  I  )  -  (virha,  q 

o  1 


2  2J 

) 


vh  (a^qir) 

2~ 

w  q  f 
o 


2  2  2 


+  4 


,  virh  a  j  q  \  1 

'  3w  fq2  '  J 

o 


2  ,  -hva.q  ,  2 

997 it  -936  ! 


13000 


(174) 


The  phase  angle  is 


*  = 


ti  +"*n 


-  -1  ( 

=  sin 


virh  a  2q2 

~  T~ 

3  w  q  l 
o 


4v  h  a  j  q 
3  i 
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(175) 
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1757 

2000 


VTrha  2r,2 
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2 


3  w 

o 


2 


VTrh  a 


1 


2  2 
"n 


2 


Again  we  see  that  if  v  =  0,  equation  (174)  is  reduced  to 
equation  (98),  i.e.,  the  response  equation  of  the  undamped  case, 
but  due  to  the  addition  of  damping,  the  response  curves  are  closed 
at  top  now  (figures  14,  15,  16). 

Up  to  the  second  approximation,  the  solutions  of  w  and  u  are  : 


w 


.  ,  ,  ,  .  2  sin  3tt£ 

-  a^qsinTT^  cos  (t-iJj)  +  w  rj  — ^  COS  (  T 


+  > 


3vhn  2  2  .  3 

'  Yoo7  ai  11  Sin  3tT^  (  4  cos  T  +  sin  T> 

vhir  2  2  . 

-  a  j  q  cos  2^4  sin  3t  (176) 

and 

—  1  2  1 

u  =  -  g  (a  piO  (4  +  —  sin  Z*£)  ( 1  +  cos  2t)  (177) 


Again  it  is  seen  that  due  to  the  effect  of  damping,  the  expression 
of  w  is  changed.  The  change  is  in  the  form  of  the  phase  angle  ^  and 
some  additional  terms  with  the  damping  parameter  v.  The  expression 
for  u,  however,  remains  the  same. 
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SECTION  6 


NUMERICAL  RESULTS  AND 
EXPERIMENTAL  VERIFICATION, 

Calculations  for  two  sample  beams  are  made.  One  is  of  the 


very  slender  type 


A! 


=  1110  ,  the  other  is  of  the 


moderately  slender  type 


A! 


r  |  I 

are  shown  in  figure  10  through  figure  17. 


=  69.3 


Response  curves 


Steel  beams  are  used.  The  physical  constants  for  steel  are: 
E  =  30  x  10^  psi,  G  =  12  x  10^  psi,  density  T  =  0.283  lb/ in.  For 
a  rectangular  beam,  the  shear  distortion  constant  k1  =  0.833,  the 
damping  parameter  v  =  18.6  ^3oj  .  Except  the  case  of  figure  17, 
all  the  other  cases  are  calculated  for  the  external  excitation  in  the 
form  of  support  accelerations. 

The  dimensions  are: 

Beam  (a)  --  very  slender  type 


length 

depth 

width 


i  =  10  inch 

h  =  — inch 
32 


b  =  —  inch 

2 

I  -6 

—  =  0.813  x  10 


Af 


The  first  natural  linear  frequency  is 


2ir 


IT 

2 


gEl 

4 

r  Ai 


=  28.5  cps 


77 


Note:  Although  the  dimensions  indicate  a  strip  type  plate 
rather  than  a  beam,  tests  show  that  the  natural 
frequency  is  equal  to  the  value  (28.5  cps )  calculated 
from  the  beam  formula. 

Beam  {b )  --  moderately  slender  type 

t  =  10  inch 

h  =  inch 

b  =  ■—  inch 
2 

I  -4 

-=—  =  2.085  x  10 

AiZ 

“l 

—  =  458.5  cps 
2it 

Note  that  the  values  of  — —  for  both  beams  are  small.  The 

AT 

correction  due  to  shear  deformation  and  rotatory  inertia  can  be 
neglected . 

Figure  10  shows  response  curves  of  Beam  (a)  for  the 
immovable  hinged  ends  case.  The  hard  spring  effect  is  very 
strong  (equation  (7  5)  ). 

Figure  1  1  shows  the  response  curve  of  Beam  (b)  with 
immovable  hinged  ends.  The  peak  amplitude  is  twice  as  high  as 
the  depth  of  the  beam.  The  hard  spring  effect  due  to  large  amplitude 
is  so  overwhelmingly  strong  that  the  influence  of  damping  is 
relatively  insignificant.  However,  the  damping  does  make  the  top 
closed.  In  calculating  the  response  curves,  equation  (7  5)  (undamped) 
and  equations  (145)  and  (158)  (damped)  are  used. 

Figure  12  shows  the  response  of  Beam  (b)  with  immovable 
hinged  ends.  But  this  time  the  peak  amplitude  is  limited  to  1/8  inch 
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{depth  of  beam  =  1/2  inch).  Equation  (75)  (undamped)  and 
equation  (158)  (damped)  are  used. 

Figure  13  shows  the  response  of  Beam  (b)  while  simply 
supported.  The  soft  spring  effect  can  be  seen  (on  the  stem  curve) 
but  not  significant. 

For  the  response  curves  of  figures  13,  14,  15  and  16, 
equations  (98)  (undamped)  and  (174)  (damped)  are  used 

If  instead  of  the  support  acceleration,  the  beam  is  excited  by 
a  periodic  uniform  loading,  the  right  side  of  equation  (56)  is  zero, 
and  consequently  the  last  term  of  the  response  equation  (98)  is 
dropped.  The  difference  is  shown  up  in  figure  17,  where  by 
equivalent  uniform  loading  means  that  if  only  the  linear  theory  is 
considered,  both  loadings  would  give  exactly  the  same  response 


relationship. 

Tests  have  been  conducted  on  the  simply  supported  very 
slender  beam.  Both  the  soft  spring  effect  and  jump  phenomena 
are  demonstrated.  Some  of  the  test  results  are  plotted  in 
figures  14,  15,  and  16- 

The  specific  damping  of  a  material  can  be  expressed  by 
D  =  J  (rm  (in-lb  per  cubic  inch  per  cycle) 
where  <r  is  the  amplitude  of  the  stress.  But  D  is  equal  to  the  area 

of  the  hysteresis  loop,  i.e., 


D  =/e 


(?-Ir)d€  (Figure  6) 


By  the  use  of  equation  (104),  the  integration  gives 

D  =  |  vECp  <C0=  WE) 


For  any  value  off  ,  v  is  determined  by  equating  these  two  expressions 

°  *  - 
of  D.  For  example,  for  Sandvik  steel  we  have:  m  =  2,3  and  J-8.9  x  10 

-3 

If  £  =  3,08  x  10  ,  v  =  2.03.  For  comparison,  response  curves  for 

v  =  2.03  and  v  =  18.6  are  both  shown  in  Figure  16. 


*See  ^26^j  /  equation  36.9,  page  36-7  and  Table  36.5,  page  36-35. 
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SECTION  7 


CONCLUSIONS 


(a)  The  nonlinear  effects  due  to  large  amplitude  on  the  response 
depends  on  support  conditions.  With  immovable  hinged  ends 
the  response  is  of  hard  spring  type;  and  for  the  case  of 

s imply ** supported  ends  the  response  is  of  the  soft  spring  type. 

(b)  Nonlinear  material  damping  in  addition  to  limiting  the  amplitude 
at  resonance  to  a  finite  value ,  modifies  the  shape  of  the 
response  curve.  It  makes  those  response  curves  of  the  soft 
spring  type  more  soft;  and  those  of  the  hard  spring  type  less 
hard. 


(c)  The  effect  of  shear  distortion  and  rotatory  inertia  on  the 


linear  or  nonlinear  response  curves  is  related  to  the  slenderness 

I  ....  I 

parameter  — - 


For  beams  with  small 


values,  the 


A!  Al 

correction  due  to  shear  distortion  and  rotatory  inertia  is 


usually  negligible. 
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Figure  10  -  Response  Curve  -  -  Very  Slender  Beam  with  Immovable  Hinged  Ends,  Undamped  Case  Only 


Loading:  w  =  5x10  inch 
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Figure  11  -  Response  Curve  --  Moderately  Slender  Beam  with  Immovable  Hinged  End 
with  Amplitude  Large  Compared  with  Depth  of  the  Beam 
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Figure  1Z  -  Response  Curve  --  Moderately  Slender  Beam 
with  Immovable  Hinged  Ends,  with  Amplitude  Comparable 
to  the  Depth  of  the  Beam 


Figure  13  -  Response  Curve  --  Moderately  Slender 
Simply  Supported  Beam 
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Amplitude  _  Amplitude 


Figure  14  -  Response  Curve  --  Very  Slender  Simply 
Supported  Beam  (with  Experimental  Results) 
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Figure  15  -  Response  Curve  --  Very  Slender  Simply 
Supported  Beam  (with  Experimental  Results) 
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igure  16  Response  Curve  --  Very  Slender  Simply  Supported  Beam  (with  Experimental  Results) 


Shaker  Loading 
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Figure  17  -  Comparison  of  Response  Curve  of  a  Very  Slender  Simply  Supported  Beam 
Due  to  (a)  Shaker  Loading,  fb)  Equivalent  Uniform  Loading 
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